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Abstract. We severely criticize the consuetudinary anal- 
ysis of the afterglows of gamma-ray bursts (GRBs) in the 
conical-ejection fireball scenarios. We argue that, instead, 
recent observations imply that the long-duration GRBs 
and their afterglows are produced by highly relativistic 
jets of cannonballs (CBs) emitted in supernova explosions. 
The CBs are heated by their collision with the supernova 
shell. The GRB is the boosted surface radiation the CBs 
emit as they reach the transparent outskirts of the shell. 
The exiting CBs further decelerate by sweeping up inter- 
stellar matter (ISM). The early X-ray afterglow is domi- 
nated by thermal bremsstrahlung from the cooling CBs, 
the optical afterglow by synchrotron radiation from the 
ISM electrons swept up by the CBs. We show that this 
model fits simply and remarkably well all the measured 
optical afterglows of the 15 GRBs with known redshift, in- 
cluding that of GRB 990123, for which unusually prompt 
data are available. We demonstrate that GRB 980425 was 
a normal GRB produced by SN1998bw, with standard X- 
ray and optical afterglows. We find that the very pecu- 
liar afterglow of GRB 970508 can be explained if its CBs 
encountered a significant jump in density as they moved 
through the ISM. The afterglows of the nearest 8 of the 
known-redshift GRBs show various degrees of evidence for 
an association with a supernova akin to SN1998bw. In all 
other cases such an association, even if present, would have 
been undetectable with the best current photometric sen- 
sitivities. This gives strong support to the proposition that 
most, maybe all, of the long-duration GRBs are associated 
with supernovae. Though our emphasis is on optical af- 
terglows, we also provide an excellent description of X-ray 
afterglows. 



Key words: gamma rays bursts, supernovae, optical af- 
terglow. X-ray afterglow 

1. Introduction 

Our information about the once totally mysterious 
gamma-ray bursts (GRBs) has increased spectacularly in 
the past few years. The rapid directional localization of 



gamma-ray bursts by the satellites BeppoSAX (e.g. Costa 
et al. 1997), Rossi (e.g. Levine et al. 1996) and by the 
Inter-Planetary Network (IPN) of the spacecrafts Ulysses, 
Konus-Wind and NEAR (e.g. CHne et al. 1999) led to a 
flurry of progress: the discovery of long-duration GRB af- 
terglows (Costa et al. 1997; Van Paradijs et al. 1998); the 
discovery of the GRBs' host galaxies (Sahu et al. 1997a); 
the measurement of their redshifts (Metzger et al. 1997b) 
that verified their cosmological origin (e.g. Paczynski 
1986; Meegan et al. 1992); the identification of their birth- 
places — mainly star formation regions in normal galaxies 
(e.g. Paczynski 1998; Holland and Hjorth 1999)— and the 
first evidence for a possible association between GRBs and 
supernova explosions (Galama et al. 1998a). 

The enormous isotropic energies inferred from the red- 
shifts and fluences of GRBs and their short-time variabil- 
ity have indicated that the GRBs must be produced by 
gravitational stellar collapse (Goodman et al. 1987, Dar 
et al. 1992). The prevalent belief is that they are gen- 
erated by synchrotron emission from relativistic fireballs 
produced by mergers of compact stars (Paczynski 1986; 
Goodman 1986; Goodman et al. 1987), by hypernova ex- 
plosions (Paczynski 1998), or by relativistic "firecones" 
(e.g. Rhoads 1997, 1999) from coUapsars or failed super- 
novae (Woosley 1993; Woosley & MacFadyen 1999; Mac- 
Fadyen & Woosley 1999; MacFadyen et al. 2001). But var- 
ious observations, including in particular the ones we shall 
extensively discuss here, strongly suggest that most of 
the long-duration GRBs are produced in supernova events 
(Dar & De Riijula 2000a and references therein) by highly 
collimated superluminal jets (e.g. Shaviv & Dar 1995; Dar 
1997; Dar & Plaga 1999). 

Various authors (e.g. Rhoads 1997, 1999; MacFadyen 
and Woosley 1999; Sari et al. 1999) have merged the no- 
tion that GRBs are produced by highly relativistic jets 
(e.g. Brainerd 1992; Woosley 1993; Shaviv & Dar 1995; 
Dar 1997; Dar 1998a; Dar & Plaga 1999) with the popu- 
lar fireball models of GRBs (see, e.g. Piran 1999 and ref- 
erences therein) to morphe the concept of "firecones" or 
similar denominations. Firecone considerations are used to 
analyze "breaks" in GRB afterglows (Rhoads 1997, 1999; 
Sari et al. 1999; Frail et al. 2001) and to extract properties 
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of the GRB engine and ejecta (Frail ct al. 2001, and ref- 
erences tlierein). The fireball idea that all radiation from 
GRBs originates from colliding shocks is certainly inter- 
esting and worth studying. But, concerning the evolution 
of AGs, the idea remains essentially untested, given the 
cavalier treatment it has received in much of the recent 
literature. In Section 2 we explain this harsh opinion. 

In recent papers the idea of jetted, supernova- 
associated GRBs was made entirely explicit with the in- 
troduction of a relativistic cannonball (CB) model of GRB 
production in supernova explosions (Dar & De Riijula 
2000a, 2000b, 2001a, 2001b, hereafter DD2000a, etc.). The 
CB model is completely different from the firecone sce- 
narios, as we explain in Section 2. The CB model, we con- 
tend, explains the main observed features of long-duration 
GRBs and of their afterglows. In particular, in DD2000b, 
we have demonstrated that the CB model predicts the 
temporal and spectral properties of the bursts of 7-rays 
correctly. 

In this paper we derive the detailed predictions of the 
CB model for the GRB optical afterglows (AGs), which 
we only sketched in DD2000a. We compare the predic- 
tions, which are analytic in fair approximations, with the 
observed optical AG of all the GRBs with known redshift. 
We show that the CB model describes remarkably well 
these optical AGs, as well as the measured X-ray AGs of 
these GRBs. 

Our detailed analysis of the AGs allows us to show 
how the nearest eight GRBs with measured redshift show 
varying degrees of evidence of an association with a su- 
pernova (SN) akin to SN1998bw: superimposed on the 
smooth AG of these GRBs one can discern the light curve 
of SN1998bw, adequately translated and red-shifted in lu- 
minosity distance, time-dependence and spectrum (e.g. 
Dar 1998b; DD2000a). In all other cases, either there are 
no late-time measurements of the optical AG, or the SN 
contribution is too dim to be resolved from the late GRB 
afterglow or the host galaxy light, even by the HST or the 
most powerful ground-based telescopes. 

In spite of the fact that we use similar vocabulary in 
what concerns the GRB engine (which is quite irrelevant 
to the AG properties discussed here) the CB model is 
completely different from the coUapsar model of GRBs 
(Woosley 1993; Woosley & MacFadyen 1999; MacFadyen 
and Woosley 1999; MacFadyen et al. 2001). The crucial 
differences are explained in DD2000a and DD2000b. One 
of them is that we use empirical facts as our guiding line, 
rather than the results of simulations that fail to explode 
SNe and lack a proper treatment of features that are no 
doubt relevant (relativistic dynamics, angular momentum 
and its transport, magnetic fields). Concerning AGs, the 
CB model does have a predecessor, the "plasmoid model" , 
of Chiang & Dermer (1997). The CB model, however, dif- 
fers in many crucial details; it is more complete and, unlike 
the plasmoid model, it is very successful in describing the 
observations. 



We devote sections 3 to 6, to make the paper self- 
contained, to a brief review of the CB model. The novel 
theoretical core of the paper is in sections 7 to 10. We 
propose there, in particular, a simple mechanism govern- 
ing the pace of radial expansion of a CB. The mecha- 
nism naturally explains the fact, observed in quasar and 
microquasar ejections, that cannonballs, faithful to their 
name, essentially stop expanding at some point of their 
voyage. We also derive predictions, which turn out to be 
in disagreement with observations, for the case of CBs that 
would continuously expand in an inertial manner. The re- 
maining sections are devoted to the description of the AG 
data in the CB model, to a detailed comparison with ob- 
servations, which turns out to be extremely successful, and 
to the extraction of conclusions. 

Three GRBs deserve special mention in their CB 
model interpretation. GRB 980425 turns out to be en- 
tirely "normal" , it is uncommonly near (z = 0.0085) but 
its emitted CBs are observed at an unusually large an- 
gle, giving it a normal 7-ray fluence. This interpretation 
(DD2000a) is strengthened by the fact that we success- 
fully predict its optical AG (dominated up to 20 months 
by SN1998bw) from its X-ray AG, which is entirely due to 
the CBs, and normal. GRB 970508 has an extremely pe- 
culiar AG light curve, which can be easily explained, but 
only if its CBs encountered a significant jump in density 
as they moved through the ISM. There are uncommonly 
early data on the optical AG of GRB 990123, which fit the 
expectations for CBs that are moving through the wind 
of the parent star, in its Wolf-Rayet phase. 

2. Uses and abuses of fireballs 

In this section we place the CB model of GRB afterglows 
in the perspective of the generally accepted views on the 
subject, which are based on the "fireball" model and on 
modifications thereof. This serves a double purpose: it 
clarifies how completely different the CB model is from 
the fireball ones, and it shows how unconvincing a frac- 
tion of the fireball literature is. 

In the fireball model, reviewed in Piran (1999, 2000) 
and Meszaros (2001), both the 7-rays and the AG of a 
GRB are made by synchrotron radiation in inward- and 
outward-moving shocks, which are produced as relativis- 
tically expanding shells collide with each other and with 
interstellar material. The attitude is often espoused that 
the actual engine producing these colliding ejecta need 
not always be explicitly discussed. The possibility that the 
ejecta may not be spherically distributed has been repeat- 
edly studied in the literature, but, in the fireball model, 
this was not done in detail prior to the influential papers 
of Rhoads (1997, 1999), who predicted abrupt breaks in 
the power law of the AG light curves. With the advent of 
GRB 990123, with its record equivalent spherical energy 
(see Table I) and an AG light curve through which it is 
possible to draw a broken power law (e.g. Figs. 1-4 of Hoi- 



S. Dado, A. Dar & A. De Riijula: Afterglows of GRBs 



land et al. 2000a) the fireball advocates (see, e.g. Frail et 
al. 2001) adopted the arguments in favour of coUimated 
GRBs (e.g. in the case of GRBs from quasars, Brainerd 
1992; in the case of a funnel in an explosion, Meszaros & 
Rees 1992; and in the case of jets in gravitational collapses, 
Shaviv & Dar 1995, Dar 1997, Dar 1998a; Dar & Plaga 
1999, DD2000a and references therein). So have fireballs 
evolved into "coUimated fireballs" , "firecones" or "conical 
fireball jets" , while maintaining the "fire" lineage. 

Consider first a proper (i.e. spherical) fireball expand- 
ing in a homogeneous (or spherically symmetric) medium. 
A conical section of this fireball would expand as a fixed- 
angle cone: a firecone. Consider next material that is 
ejected with a conical distribution. If the cone expands 
laterally at a transverse velocity vt, its opening angle, as 
viewed from the origin of the ejecta, increases with time. 
As illustrated in Fig. (Il|a) the edges of the material de- 
scribe a trumpet-shaped curve, not a fixed-aperture cone, 
as some of the names given to it may induce one to think. 
We call these "firecones" or "conical fireball jets" firetrum- 
pets, since, for vt =/= 0, that is what they are. 

Let 7(t) be the Lorentz factor of the ejecta, that di- 
minishes with time as they collide with ambient material. 
The light emitted by an element of a firetrumpet's surface 
is coUimated by its motion into an angle of aperture I/7. 
If vt/7 were constant, the firetrumpet's opening angle 
would vary as: 



9.(t)==^j(0) + ^. 



(1) 



At the time t = tb at which 6'j(t) = l/7(t), the angle of 
emission of light becomes broader than the angle of the 
cone (Rhoads 1997, 1999). Thereafter the forward light- 
collimation is less efficient, an on-axis observer would see 
up to the edge of the cone, and no longer an increasing 
fraction of the ejecta (Meszaros et al. 1999). At early times 
t < toxp, the lateral expansion of the firetrumpet may not 
be important and the ejecta's deceleration as it plunges 
through constant-density material results — as it would for 
a fixed-angle cone — in 7(r) ex r~'^/^, with r(t) the trav- 
elled distance, while at later times 7(r) oc exp[— r/roxp]- 
Rhoads assumes that these two transitions occur at the 
same time (tb = tcxp) and that they are abrupt, lead- 
ing to a break: a sudden increase in the index a of an 
AGs' power-law evolution, F oc t^". The break-time is es- 
timated as the time at which a cone which is not laterally 
expanding decelerates to 7(t) = I/60. 

Sari et al. (1999) change some of the parameters used 
by Rhoads (notably vt = c/vo to vt = c) and invert 
tb(6'j(0)) to obtain 6*3 (tb): 



0j(tb)^0j(O) = O.l 



tb ni 



6.2 h E52 



(2) 



While these theoretical developments were taking 
place, more than a dozen (mainly R-band) optical AGs 
were being observed. They did not have abrupt breaks. 
The observers (GRB 990123: Fruchter et al. 1999a, Castro- 
Tirado et al. 1999b, Kulkarni et al. 1999; GRB 990510: 
Stanek et al. 1999, Harrison et al. 1999, Israel et al. 1999; 
GRB 990705: Masetti et al. 2000; GRB 991208: Castro- 
Tirado et al. 2001; GRB 991216: Halpern et al. 2000a; 
GRB 000301c: Sagar et al. 2000b, Jensen et al. 2000; 
GRB 000418: Berger et al. 2001; GRB 000926: Fynbo 
et al. 2001, Sagar et al. 2001a; GRB 010222: Masetti et 
al. 2001, Stanek et al. 2001) fitted the slow steepening of 
AG fiucnccs to phenomenological formulae such as: 



F. = 



2F 



[(t/tb 



(t/tt 



(3) 



where ni is the local density in cm ^ and E52 is the ejecta's 
energy in 10^^ erg units. 



which interpolate between two power laws with a tun- 
able "abruptness" s, often set at s = 1. The values 
of tb extracted from these fits, and their distributions, 
have no clear meaning, since different groups use differ- 
ent parametrizations, and none of them is theoretically 
justified. 

Moderski et al. (2000), Huang at al. (2000a,b), Kumar 
& Panaitescu ( 2000) and Panaitescu & Kumar (2001) 
have modelled the light emitted by a firetrumpet without 
some of the approximations introduced by Rhoads. The 
evolution of the ejecta's deceleration is treated continu- 
ously. The emission is computed from isochronous points 
in the firetrumpet, so that light simultaneously received 
is light that had been simultaneously emitted (lifting the 
prior "approximation" that the speed of light is infinite). 
Not having an abrupt break put in by hand, no abrupt 
break is predicted. The fair conclusion is that the light 
curves are too smooth to allow for a determination of a 
break time tb (Moderski et al. 2000). 

All the firetrumpet advocates place the observer pre- 
cisely on the jet's axis, as in Fig. (yja), for no stated reason. 
It is obvious that the viewing angle is a relevant parameter 
that cannot be unceremoniously dismissed. In particular 
a non-vanishing viewing angle would contribute to erase 
even further any trace of a sharp break. Moreover, a distri- 
bution of viewing angles would completely erase a possible 
meaning to the distribution of specific tb values extracted 
from expressions such as Eqs. (g) and (||). 

The hypothetical firetrumpet ejecta behave in a differ- 
ent way from most of the highly relativistic jets observed 
in quasars (e.g. radio jets: Bridle 2000; optical jets: Crane 
et al. 1993; X-ray jets: Wilson et al. 2000) and micro- 
quasars (e.g. Mirabel & Rodriguez 1994, 1999). The ejecta 
of the real jets, as seen from their emission point up to 
the point where they eventually stop and expand, gener- 
ally subtend angles that decrease with time, exactly the 
opposite of the assumed firetrumpet behaviour of Eq. (|l]) 
and Fig. (|^a). In the analysis of these real objects (e.g. 
Pearsons & Zensus 1987; Mirabel & Rodriguez 1994, 1999; 
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Ghisellini & Celotti 2001) it is the angle of observation — 
and not the angle subtended by the ejecta — that plays a 
key role. 

Frail et al. (2001) use the published values of tb — fit 
to expressions such as Eq. (||) — to extract a set of values 
of 6'j(0). In so doing they use a modified Eq. (g), in which 
the dependence on redshift and on the efficiency of light 
production are not overlooked. In this way they reach a 
series of conclusions that their analysis does not justify. 

The firetrumpet model may, to some extent, be correct. 
We have seen that, alas, the consequences of its basic as- 
sumptions have not been properly extracted. In particular, 
the "anthropo-axial" view that the ejecta of the observed 
AGs always point to the observer has not been shown to 
be a fair approximation. 

In Fig. (|l|b) we illustrate the geometry of the CB 
model. In the AG phase, the CBs are expanding very 
slowly, or not at all (DD2000a), like the observed ejecta 
in quasars and microquasars. In contradistinction to the 
firetrumpet case of Eq. (|^), the angle with which the CBs 
are viewed from the origin decreases with time. But this 
angle is irrelevant, and negligible relative to the opening 
angle I/7 of the emitted light. The angle at which the 
ejected CBs are viewed is obviously relevant and we do 
not set it to zero by fiat. 



3. The cannonball model of GRBs 

In the CB model, long-duration GRBs and their AGs are 
produced in core collapse supernovae by jets of highly rel- 
ativistic "cannonballs" that pierce through the supernova 
shell. The detailed model is based essentially on the fol- 
lowing analogies, hypotheses and explicit calculations: 



3.1. Relativistic jets in astrophysics 

Astrophysical systems, such as quasars and microquasars, 
in which periods of intense accretion onto a compact mas- 
sive object occur, emit highly collimated relativistic jets 
of plasma. The Lorentz factor 7 = l/-y/l — v^/c^ of these 
jets ranges from mildly relativistic: 7 ^ 2.55 for GRS 
1915-^105 (Mirabel & Rodriguez 1994, 1999), to quite rel- 
ativistic: 7 = (10) for typical quasars (e.g. Ghisellini et 
al. 1993), and even to highly relativistic: 7 ~ 10'^ for PKS 
0405-385 (Kedziora-Chudczer et al. 1997). These jets are 
not continuous streams of matter, but consist of individ- 
ual blobs, or "cannonballs". (e.g. Kraft et al. 2001). The 
mechanism producing these surprisingly energetic and col- 
limated emissions is not understood, but it seems to op- 
erate pervasively in nature (In Section 7 we propose a 
mechanism capable of collimating CBs). We assume the 
CBs to be composed of ordinary "baryonic" matter (as 
opposed to e"*" e~ pairs) , as is the case in the microquasar 
SS 433, from which Lyo, and metal K^ lines have been 
detected (e.g. Margon 1984, Kotani et al. 1996). 



3.2. The GRB/SN association 

There is mounting evidence for an association of super- 
nova explosions of type Ib/Ic and GRBs (DD2000a). The 
first example was GRB 980425 (Sofhtta et al. 1998; Kip- 
pen et al. 1998), within whose error circle SN1998bw was 
soon detected optically (Galama et al. 1998a) and at radio 
frequencies (Kulkarni et al. 1998a). The chance probabil- 
ity for a spatial and temporal coincidence is less than 10~^ 
(e.g. Galama et al. 1998a), or much smaller if the revised 
BeppoSAX position (e.g. Plan et al. 1999) is used in the 
estimate (we shall show that the observed X-rays origi- 
nate in the CBs of this GRB, and not on the associated 
SN). The unusual radio (Kulkarni et al. 1998a; Wieringa 
et al. 1999) and optical (Galama et al. 1998a; Iwamoto 
et al. 1998) properties of SN1998bw support this associa- 
tion. The exceptionally small fluence and redshift of GRB 
980425 make this event very peculiar, though not in the 
CB model (DD2000a). The energy supply in a SN event 
similar to SN 1998bw is too small to accommodate the flu- 
ence of cosmological GRBs, unless their 7-rays are highly 
beamed. SN 1998bw is a pecuHar supernova, perhaps be- 
cause it is observed close to the axis of its GRB emission. 
Evidence for a SN1998bw-like contribution to a GRB 
afterglow (Dar 1999a; Castro-Tirado & Gorosabcl 1999) 
was first found by Bloom et al. (1999a) for GRB 980326, 
but the unknown redshift prevented a quantitative analy- 
sis. The AG of GRB 970228 (located at redshift z = 0.695) 
appears to be overtaken by a light curve akin to that 
of SN1998bw (located at Zbw = 0.0085), when properly 
scaled by their differing redshifts (Dar 1999b; Reichart 
1999; Galama et al. 2000). Evidence of similar associa- 
tions was found for GRB 990712 (Hjorth ct al. 2000a; 
Sahu ct al. 2000; Bjornsson ct al. 2001), GRB 980703 (Hol- 
land et al. 2000b), GRB 000418 (DD2000a), GRB 991208 
(Castro-Tirado et al. 2001) and GRB 990510 (Sokolov et 
al. 2001a). For the remaining cases in Table I, correspond- 
ing to GRBs with larger redshifts, either no late observa- 
tions of the AG are available, or the expected SN bump 
is an unobservably small effect. These conclusions will be 
strengthened by our detailed analysis of AGs in the CB 
model. All of the nearest GRBs with measured redshifts 
show various degrees of evidence for a supernova in their 
AG, suggesting the possibility of an association of all the 
long-duration GRBs with core-collapse SNe. 

3.3. The SN/GRB association 

By a SN/GRB association — as opposed to the GRB/SN 
association we have discussed — we mean the converse 
statement to that ending the last subsection: that most 
SNe of certain relatively frequent types may be associated 
with GRBs. This appears at first sight to be entirely un- 
tenable. The total rate of type Il/Ib/Ic SNe has been esti- 
mated from their observed rate in the local Universe (e.g. 
Van den Bergh & Tammann 1991) and the star formation 
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rate as function of redshift, to be Rsn — 12 ± 5s^^, or 
RsN ~ 3.8 X 10* y'^, in the observable Universe (Madau 
et al. 1998), while the observed rate of GRBs is a mere 
R-GRB — 10^ y~^ (see, for example, Lamb 2001). 

The bolometric energy fluence from a CB moving with 
a Lorentz factor 7 ::^ 1 and seen by a stationary observer 
at an angle 9 <^ 1 relative to the CB's direction of motion 
(e.g. DD 2000a) is: 

(4) 



dF 



27 



1 + 7 



2fl2 



Barring the case of GRB 980425 (whose exceptionality 
and interpretation we shall discuss) the equivalent spheri- 
cal energies of the GRBs with measured redshifts, as listed 
in Table I, range between, approximately, 2 x 10^'' erg 
(GRB 990123) and 6.6 x lO^^^ erg (GRB 970508). The 9 
dependence is the steepest parameter dependence of the 
CB model (DD2000b). It is therefore reasonable to at- 
tribute the range of observed equivalent spherical energies 
mainly to the 9 dependence (as if GRBs were otherwise 
approximately standard candles) . The observed spread in 
equivalent energy then corresponds, according to Eq. (0), 
to a spread of viewing angles between 9^0 and 9 w 2.4/7. 
Thus the geometrical fraction of GRBs which are observ- 
able (with the current or past sensitivity) is approximately 
f(7) = 27r0^/(47r) « 2.84/7^, where we have taken two 
jets of CBs per event. Compare Rsn and Rqrb to con- 
clude that an approximately one-to-one GRB/SN associ- 
ation would require beaming into a solid angle that is a 
fraction ^2.8 x 10^^ of Aw. For CBs moving with 7^ 10^, 
1(7) — 2.84 X 10~^: precisely the required beaming factor. 
That is, for a one-to-one SN/GRB association: 



10^ 



y 
7 

Moreover, 



(5) 



if the recent 



Rgrb = f(7) Rsn - (1082 ± 450) 

in agreement with observation, 
claims that the ~ (l+z)"^ dependence of the star formation 
rate continues to z > 1 (Fenimore & Ramirez- Ruiz 2000; 
Ramirez-Ruiz and Fenimore 2000; Reichart et al. 2001) 
were correct, SN of types Ib/Ic would by themselves suf- 
fice to explain the observed GRB rate. Thus, relativistic 
beaming solves the energy crisis of GRBs and may allow 
an approximately one-to-one SNIc/GRB association (Dar 
1999b; Dar & Plaga 1999; DD2000a). The above consid- 
erations leading to a GRB/SN association that may be 
as biunivocal as indicated by Eq. (^ are weakened by the 
fact that we have not taken into account effects such as the 
efficiency of GRB identification as a function of fiuence. 
It is clear, however, that for the high beaming factors we 
have advocated, the GRB and SN rates are quite com- 
parable. In the CB model, the "special" character of SN 
1998bw is due to the fact that it is observed very close to 
the GRB axis. 

Previous analyses of GRB/SN associations, except 
that in DD2000a, were based on a power-law extrapola- 
tion to late times of the early-time GRB afterglows. Here 



we shall use the CB model, instead, to calculate the GRB- 
afterglow light curves at all times. This procedure leads, as 
we shall see, to a much better exposition of the GRB/SN 
association. 



3.4- The GRB engine 

We assume that in core-collapse SN events a fraction of the 
parent star's material, external to the newly-born compact 
object, falls back in a time very roughly of the order of one 
day (De Riijula 1987, DD2000a). Given the considerable 
specific angular momentum of stars, it settles into an ac- 
cretion disk and/or torus around the compact object. The 
subsequent sudden episodes of accretion — occurring with 
a time sequence that we cannot predict — result in the 
emission of CBs. These emissions last till the reservoir of 
accreting matter is exhausted. The emitted CBs initially 
expand in the SN rest system at a speed /?c/7, with /3 c 
presumably of the same order as the speed of sound in a 
relativistic plasma (/? = l/\/3), or smaller. The solid an- 
gle a CB subtends is so extremely small that presumably 
successive CBs do not hit the same point of the outgoing 
SN shell, as they catch up with it. These considerations 
are illustrated in Fig. (g). 

3.5. The GRB 

From this point onwards, the CB model is not based on 
analogies or assumptions, but on processes whose outcome 
can be approximately worked out in an explicit manner. 
The violent collision of the CB with the SN shell heats the 
CB (which is not transparent at this point to 7's from tt'^ 
decays) to a surface temperature that, by the time the CB 
reaches the transparent outskirts of the SN shell, is ~ 150 
eV, further decreasing as the CB travels (DD2000b). The 
resulting CB surface radiation, Doppler-shifted in energy 
and forward-coUimated by the CB's fast motion, gives rise 
to an individual pulse in a GRB (DD2000b). The GRB 
light curve is an ensemble of such pulses, often overlap- 
ping one another. The energies of the individual GRB 7- 
rays, as well as their typical total fiuences, indicate CB 
Lorentz factors of 0(10^), as the SN/GRB association 
does (DD2000a). In the CB model, unlike in the shocked- 
fireball models, the photons of the GRB proper are not 
made by synchrotron radiation which, as we shall see, is 
subdominant at this stage of the evolution of a CB. 

4. Afterglow components 

In the CB model, the persistent radiation in the direction 
of an observed GRB has three origins: the ejected CBs, 
the concomitant SN explosion, and the host galaxy. These 
components are usually unresolved in the measured "GRB 
afterglows", so that the corresponding light curves and 
spectra refer to a cumulative energy flux density: 



Fag 



'CBs 



'SN 



HG : 



(6) 



S. Dado, A. Dar & A. De Riijula: Afterglows of GRBs 



with F = ;/ dN^ / (dt d:/ d A) . 

The emission of the GRB's host galaxy is usually de- 
termined from measurements at times late enough for the 
CB's afterglow and the SN light to have become com- 
paratively weak (e.g. Sokolov et al. 2001b and references 
therein). This assumes that the host galaxy's emission is 
steady on periods of a few months. There is no indication 
of GRB host-galaxy variability on such time scales. 

Core-collapse supernovae (SNII/Ib/Ic) are far from be- 
ing standard candles. But if their explosions are fairly 
asymmetric — as they would be if a fair fraction of them 
emitted jets of CBs — much of the variability could be a 
reflection of the varying angles from which we see their 
non-spherically expanding shells. Exploiting this possibil- 
ity to its extreme, we shall use SN1998bw as an ansatz 
standard candle, associated with every GRB of known z 
(Dar 1999b; DD2000a). The adequacy of this bold hypoth- 
esis can be judged from its rather surprising success. 

Let the energy flux density of SN1998bw be Fbw[t^, t]. 
For a similar SN placed at a redshift z: 



FsNht] = 



1 + z D2(zbw) 



1 + 

^ bw 



Zbw Dl(z) 
1 + z 



1 +Zbv, 



1 + Zbw 



1 



A(i/,z)^ 



(7) 



where Dl(z) is the luminosity distance^ and A(j/, z) is the 
extinction along the line of sight. The extinction in our 
Galaxy is reasonably well measured, but for the GRBs' 
environments it must be estimated from the spectra of 
each particular AG and host galaxy. 

The contribution of CBs to the GRB afterglows re- 
quires a much more detailed discussion. 

5. Times and frequencies 

Four "clocks" ticking at different paces and three different 
scales of frequency need be considered in the cannonball 
model of GRBs and their afterglows. 

Let 7 = 1/ y/l - l3^ = Ecb/(Mcbc2) be the Lorentz 
factor of a CB, which diminishes with time as the CB hits 
the SN shell and as it subsequently ploughs through the 
interstellar medium. Let tsN be the local time in the SN 
rest system, tcB the time in the CB's rest system, tob the 
time measured by a nearby observer viewing the CB at an 
angle 9 away from its direction of motion, and t the time 
measured by an earthly observer viewing the CB at the 
same angle, but from a "cosmological" distance (z ^ 0). 
Let X be the distance travelled by the CB in the SN rest 
system. The relations between the above quantities are: 



dtsN ^ 
dtcB 



dx 

7 dtcB = -3- : 
p c 

- S dtob ; 



^ The cosmological parameters we use in our calculations are: 
Ho = 65 km/(s Mpc), ^m = 0.3 and f^A = 0.7. 



dt = (l 



1 -t- z 
z) dtob = r- dtsN 



J S 



where the Doppler factor 6 is: 

1 27 



S = 



(8) 



(9) 



7(1-/3 cos 6*) (1-^02^2) ' 

and its approximate expression is valid for 9 <^ 1 and 
7 ::^ 1, the domain of interest here. Notice that for large 7 
and not large ^7, there is an enormous "relativistic aber- 
ration": dt ~ dtsN/7^, and the observer sees a long CB 
story as a film in extremely fast motion. 

The frequency of the photons radiated by a CB in its 
rest system, i^cB , their frequency in the direction 9 in the 
local SN system, z^sn, and the photon frequency v mea- 
sured by a cosmologically distant observer, are related by: 

i^CB = —r- ; vsN = {I + z) v , (10) 



with S as in Eq.(|^). 



6. The cooling of CBs 

As a CB pierces through the SN shell, its surface is heated 
by the collisions with the shell's constituents, and cools 
down from an early maximum temperature because of the 
decreasing density of the shell's material it collides with (a 
detailed description of the CB-shell collision can be found 
in DD2001b). At this early point of a CB's avatars, the 
internal radiation pressure is very large. Thus, in studying 
the properties of the 7 rays (DD2000b), we assumed the 
CBs to expand (in their rest system and at early times) 
at a speed comparable to that of sound in a relativistic 
plasma (c/\/3). This fast expansion implies that it is a 
good approximation to treat CBs, in their rest system, as 
spherical objects. 

Let Njct be the baryon or electron number of the en- 
semble of CBs in a jet, which we have estimated to be 
Njet -- 6 X 10^1 (e.g Eq. (5) of DD2001b, for 7i„ = 10"), 
which is close to that of the Earth (N^ ~ 3.6 x 10^^). On 
average, GRBs consist of 5 to 10 significant pulses, so that 
a single CB may have one order of magnitude fewer con- 
stituents. As they exit the shell and enter the interstellar 
medium (ISM) , CBs become transparent to their enclosed 
radiation when they reach a radius: 



Rt 



-— NcB ctt 

47r 



-,1/2 



(10^^ cm) 



NcB 



-,1/2 



6 X 105" 



(11) 



where <tt = 6.65 x lO^^s (,-^{2 jg (;]^g Thomson cross sec- 
tion. We can use Eqs. (||) to conclude that, if the CBs are 
expanding at a fraction /^trans of the speed of hght^, they 
reach a size Rtrans in an observer's time: 



tt. 



.CB 



(1 + Z) Rti-a 
S /3trans C 



(12) 



'^ The quantity /3trans is nearly identical to /3out, the trans- 
verse speed as the CB exits the SN shell, introduced in our 
previous work on the CB model (e.g. DD200a,b). 
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For the reference value of Nqb in Eq. (|Tl|), 7 = 1/6' = 10^ 
and /?trans = l/(3-\/3), CBs become transparent in a mere 

ttrans "^ '^■^ S. 

The GRB is emitted by the CBs from a distance of 
0(1) radiation length from the exterior of the SN shell, 
when their temperature is T^ ~ 150 eV and their ra- 
dius, for our typical parameters, is R^ ~ 2.5 x 10^^ cm 
(DD2000b, DD2001b). Soon thereafter, travelling in a thin 
environment and expanding fast, the CBs should cool in 
an approximately adiabatic way. Their temperature at 
ttrans IS then: 



Ttra 



R^ 



Rt 



■T~4.0eV. 



(13) 



^From about one third of ttrans onwards, the CBs would 
appear to be "coUisionless" to the ISM hadrons pierc- 
ing through thenfl, since the high-energy nucleon-nucleon 
cross section (ctn '^ 4 x 10~^^ cm^) is about one order of 
magnitude smaller than ax and the condition for "trans- 
parency" to the ISM particles is, up to a numerical factor 
of 0(1), analogous to Eq. (^. 



7. The expansion of CBs 

When a CB, in a matter of (observer's) seconds, becomes 
transparent to radiation, it loses its internal radiation 
pressure. If it has been expanding up to that moment at 
a speed comparable to that of relativistic sound, should it 
not inertially continue to do so? The fact that it is coUi- 
sionless makes the conclusion seem unavoidable. But the 
CBs emitted by many quasars appear, within the resolu- 
tion of the observations, not to expand laterally for most of 
their trajectory, before their forward motion nearly stops. 
What may the reason be? 

The ISM the CBs traverse has been previously par- 
tially ionized by the forward-beamed GRB radiation. The 
neutral ISM fraction is efficiently ionized by Coulomb in- 
teractions as it enters the CB. In analogy to processes 
occurring in quasar and microquasar ejections, the bulk 
of the swept-up ionized ISM particles are multiply scat- 
tered, in a "coUisionless" way, by the CBs' turbulent mag- 
netic fields. As illustrated in Fig. (|^), in the rest system 
of the CB these particles are isotropically re-emitted into 
the ISM. In the rest system of the parent SN they are for- 
ward coUimated and boosted to an energy '-^mc^7^ (Dar 
1998b). The isotropic re-emission implies an inwards force 
on the CB's surface. Assume that the bulk of the ISM par- 
ticles are not re-energizcd by the CB's turbulent fields. Let 
R be the CB's radius and let np be the proton ISM num- 
ber density. The rate at which the ISM protons impinge 
on the CB is r — 7rR^C7np, with 7np the ISM proton 
density seen from the CB's rest system. The momentum 



^ ISM particles that get entangled in the CB's magnetic field 
would not be coUisionless after such a very short time. 



(or, for large 7, the energy) of these protons isotropically 
leaving the CB is, per unit surface. 



Epc 
47rR2 



■mp7 npC 



(14) 



During the first hours after the GRB time, the CBs 
are still fully ionized and cooling rapidly by expansion 
and bremsstrahlung (DD2001a). Their full constituency of 
relatively cold ions, electrons, cosmic rays and entangled 
magnetic fields is electromagnetically coupled, and sub- 
ject to the very large inwards pressure of Eq. (|lj). This 
stabilizes the CB's radius to an asymptotic value Rmax- To 
estimate it, since the initial expansion velocity is not fully 
relativistic (/St^ans ^ 1), we may use Newton's equation: 



P = 



M, 



CB 



d^R 



47rR2 dt2,3 



(15) 



with P as in Eq. ([l4|), and integrate, to obtain, for a con- 
st antn n,,: 



R„ 



Rt 



3 NcB /3trai] 

2 7rnp72 



(16) 



where we have approximated 7 by its initial value because 
the asymptotic radius, as we shall see, is reached much 
before the CB has had the time to decelerate significantly. 
For a value np = 1 cm~^ of the ISM density close to the 
progenitor, Ncb = 6 x 10^*^, 70 = 10^, Eq. (^ gives 

R„,ax = 2.2 X 1014 cm [lO^^ cm] for Arans = 1/(3^3) 

[1/(10V3)]. 

The interval (in the CB's rest system) between the 
times when the CB has radius Rtrans and radius R can be 
deduced from Eqs. (|lj) and (|5|) to be: 



t(R) = 



dx 



'R3 



R3 

Hrans 



Rtr. 



„ C/3ti.a 



(17) 



The observer's time is shorter by a factor (1 -I- z)/S. In 
Fig. (HI) we invert Eq. (|^ to show the CB's radius as a 
function of observer's time (in minutes), for z = l,S = 10'^, 
and the other typical parameters quoted in the previous 
paragraph, for two choices of /3trans- The CB is seen to 
expand linearly at a speed close to the initial /3trans c, and 
then to settle fast into an approximately constant radius. 
To a good approximation the steady radius is reached in 
an observer's time: 



t. 



1+Z Rn 



^ Ptrans C 



(18) 



4 A density distribution falling with distance to the progen- 
itor star as x~^ (in a certain distance-domain) gives similar 
results, but in terms of more parameters. It suffices for the 
moment to use a constant value representing an average den- 
sity close to the progenitor, particularly at the considerable 
distances from the progenitor at which the value of Rmax is 
reached. 
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which yields ^ 1.2 [1.9] minutes for /3tians = l/(3-\/3) 
[l/(10-\/3)], the examples in Fig. (Q). Thus, typically, a 
few observer's minutes after the GRB, the CBs are ex- 
panding very very slowly. 

To estimate the internal magnetic field of the CB after 
it stopped expanding, Boo , we conjecture that the bulk 
of the kinetic energy of the CB's expansion after it be- 
comes transparent is converted to internal magnetic en- 
ergy. Since R^^^ » Rt^.^ns and (by hypothesis) Pl^^^ < 1, 
this means: 



iTT 3 B|^ 

3 "'"'' Stt 



NcB mp P\ 



•2 . 

P ^^trans 



which, upon substitution of Eq. dlq), yields: 



B„ 



100 



r np 1 


1/2 


r 7 1 

1 103 J 


.1 cm~3. 



Gauss 



(19) 



(20) 



The very large field of Eq. (^0|) is consistent with the fact 
that, in order to be able to sustain the inwards pressure 
of the isotropically re-emitted protons that it ejects, the 
magnetic field within the CB must have a pressure (or en- 
ergy density) comparable to the pressure P of Eq. (14). 
The condition P = B^/(87r) exactly reproduces Eq. (2C). 
We could also have added the building-up magnetic pres- 
sure to Eq. ( p^ , to obtain a result for Rmax differing from 
that of Eq. (EOh by 2-'^/3^ This affects our conclusions in- 
significantly and, in any case, we cannot pretend to have a 
detailed understanding of the magnetohydrodynamics of 
turbulent plasmas. 

We have argued that the re-emission of the ISM pro- 
tons, isotropic in the CB's rest frame, is what makes CBs 
stop expanding at a speed that ab-initio was semirela- 
tivistic. This may also explain the surprising quasar and 
microquasar observations. There may be other reasons — 
such as Coulomb-interaction ram pressure from the ambi- 
ent material — for CBs to stop expanding significantly at 
some point of their voyage. The strength of our conclusion 
that CBs expand slowly — or not at all — during the AG 
phase should be judged from the ability of the CB model 
to describe the AG observations. 



8. The dominant afterglow mechanisms 

In the CB model the GRB emission in 7 rays is mainly 
of thermal origin (although it does not have a thermal 
spectrum) and, in a fixed energy interval, it decreases ex- 
ponentially with time (DD2000b). A few seconds after the 
last GRB pulse (the last CB), this pseudothermal emission 
becomes a subdominant effect. For the next few hours, the 
evolution of a CB is interestingly complicated. In partic- 
ular, its originally ionized material should recombine into 
hydrogen and emit Lyman-a lines that are seen Doppler- 
boosted to keV energies (DD2001a). Later, the CBs settle 
down to a much simpler phase, which typically lasts for 
months, till the CBs finally stop moving relativistically. 



Because of the CBs' large Doppler factors, radio emis- 
sion in their rest frame is boosted to optical light in the ob- 
server's frame while their emitted optical light is boosted 
to the soft X-ray band. Radio emission from astrophysical 
plasmas at eV temperatures is mainly due to synchrotron 
radiation from relativistic electrons, whereas their optical 
glow is usually due to bremsstrahlung and line emission. 
For parameters in the general vicinity of the ones we have 
argued to be "typical" of the CB model, the X-ray AG is 
initially dominated by thermal bremsstrahlung (and line 
emission) and by synchrotron radiation thereafter, while 
the optical AG, generally observed later, is dominated by 
synchrotron radiation. In this section we analyse these two 
dominant mechanisms, relegating to Appendix 1 the dis- 
cussion of various subdominant ones. 

9. Thermal bremsstrahlung: the early X-ray AG 

When it becomes transparent, a CB cools down mainly 
by thermal-electron bremsstrahlung (TB) in ep collisions 
and by expansion. The comoving TB emission rate (e.g. 
Peebles 1993) is: 



Lbrc 



77 



nc^T^/^ergcm3s -^ 



(21) 



where rj = 1.435 x 10~^^, iTo is the electron density of 
the CB and, here and in the rest of this section, T is its 
temperature in Kelvin and the remaining quantities are in 
c.g.s. units. Thus, as long as the CB is fully ionized, its 
total comoving TB energy-loss rate is: 



dE, 



CB 



3N2,bT1/2 
''1 — A T,^ erg s 



dtcB ' 47rR3 

while its thermal energy is: 
EcB=^3NcBkT, 



(22) 



(23) 



with k — 1.38 X lO^^^erg/deg Boltzmann's constant. As 
long as TB dominates the cooling and for a constant ex- 
pansion rate R(tcB) ~ AransCtcB, Eqs. (|l|, |l|, H, HI) 
yield T for tcB > tg^^,: 



Ti/^(tcB) 



K 



77 N, 



CB 



K 



T. 



16 7rk/3transCR2 
1/2 V 



(24) 

12k/3transCCrT 

A distant observer at an angle 9 relative to the 
CB's direction of motion receives this radiation in a 
Doppler-boosted, coUimated and time-aberrant form (e.g. 
DD2000a). At a luminosity distance Dl(z), the total 
power (integrated over frequencies) per unit area is: 



dF 
dtdr^ 



3r;N^BT^/^^o^ 
16 7r2 R3 Dl 



erg cm s , 



(25) 



where Sq = (5[7o,0] is given by Eq. (m with the initial 7 
value, which does not change in the very short time during 
which TB dominates the X-ray AG. At the transparency 
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radius and reference baryon number of Eq. (|llL the tem- 
perature is ^ 4.0 eV ~ 4.6 x 10''^ K, as in Eq. (|13|), and for 
S = lO'^ and z = 1, the predicted energy flux, per CB, is 
1.7 X 10~* ergcm^^ s~^. This radiation's spectrum is that 
of bremsstrahlung (a flat Edn-y/dE), typically extending 
to E - 3 T 5/(1 + z) - 12 kcV, in the X-ray domain. 

For our typical CB parameters, synchrotron emission 
takes over, as we shall see in the next section, before 
the CB reaches its asymptotic radius and before the K 
term in Eq. ( p4|) becomes important. During this early 
phase (for t ~ ttrans), the TB emission by the CB's 
electrons declines with time, if cooling is dominated by 
bremsstrahlung losses, as R~^ T^/^ ~ R~^ ~ t~^. Should 
the temperature decrease be dominated by adiabatic cool- 
ing (3NcBkdT = -PdV), the substitution P = 2nckT 
results in T oc R^^ and in a TB emission declining as 
r-3ti/2^r-4_ -,--4 ^ decline as fast as t"'' or t'^ has 
been observed in early X-ray afterglow observations, e.g. 
GRB 920723: Burenin et al. 1999; GRB 970508: Piro et 
al. 1998; GRB 970828: Smith et al. 2001; GRB 990510: 
Plan et al. 2001; GRB 010222: In 't Zand et al. 2001. In 
less than a minute of observer's time, this emission mecha- 
nism is overtaken by synchrotron radiation, whose decline 
is much slower, as we proceed to discuss. 

10. Synchrotron radiation: the optical afterglow 

In this section we study the various effects resulting from 
the interaction of a CB's entangled magnetic field with the 
ISM particles that it sweeps as it travels. The CBs lose mo- 
mentum by sweeping the nuclei of the ionized ISM, and 
re-emitting them isotropically (in the CB's rest system) at 
an energy comparable to their incoming one. This allows 
us to predict the law of CB deceleration: the behaviour 
of its decreasing Lorentz factor 7(1). The incoming elec- 
trons, suffering collisions with the magnetic domains, and 
losing energy effectively by synchrotron radiation, acquire 
a predictable power-law energy spectrum, which implies 
a given distribution of the emitted photons. The emitted 
energy rate is equal to the rate at which the ISM electrons 
bring energy into the CB in its rest system; this provides 
the absolute normalization of the AG light curve. The dy- 
namical time for the energy supply by the swept-up ISM 
electrons is much longer than the time it takes the elec- 
trons to acquire a power-law energy distribution and to 
emit synchrotron radiation, justifying a quasi-steady-state 
analysis. 

The rate at which the energy of the ISM electrons en- 
ters the CB (in its rest frame) is: 



dE 



CB 



3„2 



dtcB 



TT R Uo nio c 7 



(26) 



where the incident ISM electron energy is 7 rUg c^ , and the 
extra power of 7 originates in the Lorentz contraction of 
the ISM electron density, Ug. An observer at a luminos- 
ity distance Dl and at an angle 6 relative to the CB's 



direction of motion receives a total power per unit area 
(integrated over frequencies): 



dF 
dtdfl 



ttR^ 



■^2g4 



inDl 



(27) 



At the very early time of CB transparency, for R = 
R-trans — 10^^ Cm, Uq = 1 cm^'^, and for z = 1, 
7 = (5 = 10^, the above expression yields 1.27 x 10~^ 
erg/(cm^ s), which is comparable to the bremsstrahlung 



emission of Eq. (^2|). But, as we shall see anon, the syn- 
chrotron radiation has a much softer spectrum than that 
of bremsstrahlung, and the latter mechanism dominates at 
early times in the X-ray domain. We have seen that in a 
matter of minutes the asymptotic radius Rmax of Eq. (16) 
is reached, so that thermal bremsstrahlung has decreased 
by more than three orders of magnitude, while the syn- 
chrotron radiation has increased by two or more orders of 
magnitude, to become the dominant emission mechanism 
at all frequencies. 

To estimate a "dynamical time", Tdyn, for the energy 
deposited by electrons in the CB, we forestall that the 
observed afterglow fluences are the ones expected in the 
CB model. We can then use Eqs. ( |2^ ) and ( |27| ) to deduce 
a total typical Ecb ~ 3 x lO''^ erg and to conclude: 



'^dyn 



1 dE, 



CB 



Ecb dt 



CB 



B.2 X 10^ s) 



10^ 

7 



(28) 



for R = R„iax = 2.2 x 10^"* cm and n^ = 10"-^ cm^^^ 

In the CB's rest frame, the light crossing time for the 
asymptotic radius Rmax is ^cr ~ R-max/c ^ W^ s. The time 
for electrons that enter a CB to redistribute their energy as 
they bounce off a few magnetic "sub-domains" is a fraction 
of Tcr; this is much shorter than Tdyn, so that a cosmic- 
ray-like "source" distribution of electrons (a power law 
in energy) is steadily generated. The synchrotron cooling 
time of electrons in the CB's rest system is : 



6 7rmec2 

7cCCrTB2 



10^ 

7c 



100 Gauss 



B 



(29) 



In the above equation, we have distinguished 7e (the 
Lorentz factor of an electron in the CB, in the CB's rest 
system) from 7 (the CB's bulk- motion Lorentz factor). 
Even for 7c — > 1, Tdyn ^ Tgyn. The Larmor radius of 
electrons tl — Pc/(eB) in a B = 100 Gauss magnetic 
field is rL ~ (1.5 x 10^ cm)(7e/10^), so that even for 
very high energies, the residence time of the electrons 
in the CB is much longer than Tgyn- The above inequali- 
ties imply that for electrons of all energies, a spectrum of 
"Fermi-accelerated" , radiation-loss-modulated electrons is 
steadily generated. 

The incoming ISM particles are Fermi-accelerated by 
the turbulent magnetic fields inside the CBs to a comoving 
"cosmic-ray" spectral distribution dn/dE ^ E^'^p, with 
/3p ~ 2.2, as indicated by simulations (Bendarz & Os- 
trowski 1998), analytical estimates and the interpretation 
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of the observations of cosmic rays (e.g. Dar & Plaga 1999). 
The acceleration being due to deflections by magnetic 
fields, the spectral shape of the "source" distributions of 
protons and electrons ought to be the same: /3e — /3p. The 
index of the equilibrium electron spectrum, modulated by 
radiation losses, is one unit higher: /3c — /3p + 1 ~ 3.2 
(Dar & De Riijula 2001). The emitted synchrotron radia- 
tion has a power spectrum with index a = {Pe — l)/2, so 
that: 



dF 
dv, 



dn. 



CB 

a 



di^. 



'- oc u-^ 



1.1 . 



(30) 



The spectrum of Eq. (|30|) should roughly extend between 
two cutoff frequencies, i-',nin and J^max, that reflect the en- 
ergy of the incoming electrons and of the maximally accel- 
erated ones. The integrated spectrum of Eq. ( p7| ) is pro- 
portional to t'jfjY^ and z^min oc 7c where 7c is the electrons' 
Lorentz factor above which they are in radiative equilib- 
rium. Since the individual frequencies v and the limiting 
frequency fmin all refer to the CB's rest frame, and are 
Doppler-shifted by its motion {v (x 5 i'cb) as in Eqs. (fO), 
the non-frequency-integrated version of Eq. (g^ — that 
is, the predicted spectral energy density for a GRB with 
a number n^^^ of CBs — is: 

- ~ '^^t'"' *' ^^ - n,3 (a - 1) ^ [1 + z](i-") mc c^ 



dt dv dn 

nc(x[t]) [R(t)]M7(t)]^" 
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27(t) 



3eBoc 

47r nic c 



0.42 



B. 



100 Gauss 
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(31) 



where we used the explicit form of 6, Eq. (|^), no(x[t]) is 
the density along the CB's trajectory x — J7(5cdt/(1 -|- 
z), and the overall normalization is obtained by assuming 
that relativistic electrons are in radiative equilibrium. The 
predicted normalization for a GRB is just an estimate, for 
part of the energy deposited in the CB by ISM protons, 
as well as a fraction of its magnetic energy, may also be 
emitted as synchrotron radiationn. 

When a "typical" CB, within a minute or two after the 
end of its GRB, reaches its final radius Rmax ~ 2 x 10^^ 
cm, and for Uo = 1 cni^^, z = 1, and ^ ~ 5 ~ lO'^, 
Eq. (|l|) yields 1.6 x lO'^ erg/(cm2 s), of which 6.4 x IQ-io 
erg/(cm^ s) (~3.7%) is in the 2-10 keV X-ray range and 
5.6 X 10^^ erg/(cm^ s) (~ 3.2%) is in the visible range 

^ At very low frequencies, such as those corresponding to ra- 
dio waves in the observer's frame, Eq. (pll) is expected to break 
down for a variety of reasons: a deviation of the low-energy elec- 
tron spectrum from a universal power-law, inverse synchrotron 
and inverse bremsstrhalung self-absorption, plasma frequency 
cutoff and the effect of competing mechanisms other than syn- 
chrotron radiation which all depend on the exact density profile 
and ionization state of the CB. 



(3900 A< A < 7600 A) corresponding to a spectral flux 
density of 2 Jansky (8 magnitude!) in the R band. For 
the next few hours 7(t) does not change significantly and 
the X-ray and optical AGs vary as nc(x[t]). This variation 
should in general be a decline, since the CBs are depart- 
ing from a dense region. Such a decline may have been 
observed both in the optical band in the case of GRB 
990123 (Akerlof et al. 1999), which rose above - 9th mag- 
nitude tens of seconds after the GRB's onset, not far from 
our estimate with "typical parameters" (this GRB, whose 
early optical AG we shall discuss in detail, is at z = 1.6, 
but its initial 7 and 5 are very large, see Table III). 

There may be an energy below which the synchrotron 
cooling time of the electrons in the CB is longer than 
their acceleration time. If so, the electron spectrum has 
an index /3e — /3p and the synchrotron radiation below 
a certain frequency would have an index a ~ 0.6. We 
have implicitey assumed that this frequency is below the 
smallest optically observed ones, an assumption that the 
data generally support. 

The flux of Eq. ( pi| ) depends on Vc as i'""^, roughly 
the 10-th root of it. The actual value of Vc is therefore 
quite irrelevant to the optical and X-ray AGs discussed 
here. In our study of radio AGs (Dado et al. 2002) we find 
that, in the CB model, v^ is actually the characteristic 
synchrotron frequency emitted by the electrons that enter 
the CB with a Lorentz factor (in the CB's rest frame) 
7c — 7(t), that is Vc ^ 0.22 7c :^l, with i^l oc B ex 7 the 
Larmor radius in the CB's magnetic field. The spectral 
index gradually changes from a ~ 0.5 to a ~ 1.1 at this 
frequency. For a spectrum with this transition in its power 
law, Eqs. ( pl|) are to be modified as follows: 

[7(t)]^" - [7(t)]' 

(/3p-2)(3~/3p) 
(a-1) -^ 



a « 0.5 if (1 + z)v <Siyc , 
a w 1.1 if (1 + z)iy > 6vc. 



(32) 

The difference between Eqs. (plh and (p2) is only relevant 
to the spectral shape of some very early optical AGs, and 
to the ensemble of the radio AGs (Dado et al. 2002). 

10.1. The density of the Inter- Stellar Medium 

The density of ISM protons very close to a GRB progenitor 
plays a role in determining the asymptotic radius, Rmax, 
of a CB, see Eq. (p^. The density of ISM protons along 
a CB's trajectory controls, as we shall see, the evolution 
in time of the CB's Lorentz factor 7(t). This function, 
and the density of interstellar electrons along the CB's 
trajectory, determine, via Eq. (pi]), the AG properties. 
Clearly, we must discuss these densities in some detail. 

An analysis of historical SNe in the Galaxy (e.g. Hig- 
don and Lingenfelter, 1980), of SNe in the LMC (e.g.. 
Dune et al. 2001) and of SNe in late-type galaxies (Kenni- 
cut et al. 1989; van Dyke et al. 1996; Higdon et al. 1998) 
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indicates that 85± 10% of SNe occur in superbubbles (e.g., 
Lingenfelter et al. 2001). These are spaces of typical size 
0.1 to 0.5 kpc, surrounding star-formation regions, that ex- 
tend all the way into the galactic halo, and from which the 
ISM has been swept away by massive-star winds and previ- 
ous SNe, resulting in an ISM with a low density (n ~ 10^^ 
to 10~^ cm~'^) comparable to that in a galactic halo. 

At their Wolf-Rayet phase, massive stars that finally 
produce SNell/Ib/Ic emit strong winds with typical ve- 
locities, V, close to lO^kms"""^, at a typical mass-loss rate 
M = a few 10~^ Mq yr"-'^, over the last ~ 10^ yr before the 
SN event. The density close to an imminent SN is gov- 
erned by the recent Wolf-Rayet wind and ejections, and 
declines roughly quadratically with distance as: 



M 



4 7r vx^ 
(0.18 cm"3) 



M 



lO-4M0yr-i 



lO^kms 



11 r 



1 pc 



(33) 

2 



till X'^ 10 pc, where the density becomes that of the sur- 
rounding superbubble (e.g., Ramirez- Ruiz et al. 2000 and 
references therein). 

10.2. Complications simplified 

To predict the explicit time-dependence of the AG from 
a CB, one needs to know R(t), 7(t) and the ISM density 
profile, no(x), along the CB's trajectory. Moreover, the 
various CBs that produce the different pulses in a single 
GRB have slightly different physical parameters (baryon 
number, Lorentz boost) that lead to the differences be- 
tween the individual 7-ray pulses of a given GRB. The 
large powers of the Lorentz and Doppler factors in Eq. ( pi| ) 
favour the contribution of CBs with the largest 7. Given 
the extremely small fraction of solid angle that a CB spans 
as viewed from the SN centre, we do not expect consecu- 
tive CBs to hit the SN shell on the same spot (DD2000a). 
But it is in principle possible that the initial expansion 
and slowing down of the CBs by the SN shell and the 
ISM merges several of them into a single leading CB in 
the AG phase. One seems to be faced with a plethora of 
parameters and possibilities. 

We shall find it sufficient to characterize the various 
ISM densities that the CB encounters by two constant 
densities. One of them is the average proton density very 
close to the parent star, that determines the fast-reached 
asymptotic radius of the CB. For its reference value we 
adopt Up^ = 1 cm~^. The other is the proton or electron 
density in the superbubble and in the galactic halo, for 
which we adopt as reference n = 10^"^ cm"'^. 

All other putative complications previously quoted are 
eased by the fact that the times over which AGs extend 
are much longer than the typical intervals between GRB 
pulses, so that the AG light curve is the sum of tempo- 
rally unresolved individual CB afterglows. We can there- 



fore characterize, as in Eq. (|3l|), the AG with the param- 
eters of a single CB, whose actual values would represent 
a weighted average. 

10.3. The slow down of a CB 

A CB ionizing and ploughing through an ionized ISM 
of roughly constant density, would lose momentum at a 
roughly constant rate, independent of whether the ISM 
constituents are rescattered isotropically in the CB's rest 
frame, or their mass is added to that of the CB. Energy- 
momentum conservation for a highly relativistic CB of 
initial mass Mcb — Ncb nip results in the deceleration 
law (DD2000a): 



d7 = — 



TT R2 Up 72 



dx. 



(34) 



The element dx of travelled distance (for 7^1) is, ac- 
cording to Eqs. (0), related to the observer's time interval 
dt as: 



dx = 



C7(5dt 
Tl + z) 



(35) 



The law governing the CB's expansion rate is the differ- 
ential version of Eq. (17), to wit: 



dR = C /3tra 



R3 



R3 



R3 - R3 

^max ^trans 



dt. 



(36) 



The above set of three equations can be integrated numer- 
ically for any given 70 = 7(t = 0), R(t = 0) = Rtrans and 
ISM density along the CB's path np(x). 

We limit our discussion to the case of CBs that, having 
reached their asymptotic radius, are moving through a 
constant-density medium, a case for which there are useful 
analytical expressions for the solution 7(1) of Eqs. (p4|), 
( P5| ) and (^. A constant density is a fair approximation, 
for, in a very short observer's time, the CB reaches the 
distance from the SN at which the density is that of the 
surrounding superbubble. To estimate this brief time, we 
may use the initial Lorentz and Doppler factors to obtain: 



t 



(1 + z) 
C7o(So 



X ~ (34min) 



10 pc 



(37) 



where we used the typical parameters z = 1, 70 = 10^ and 
9 = 1/70- This constant density approximation is also jus- 
tified a posteriori by the agreement between our predicted 
AG light curves and the observed ones. 

We have argued that CBs reach a steady radius Rmax 
in a few observer's minutes. To ascertain the CB's slow- 
down law for constant radius and constant ISM density, 
we may substitute Eqs. (||) and (|5|) into Eq. ( p4|) and 
integrate, to obtain: 

1 



To 


1 

.7 


1 
70. 


6ct 


(l+z)Xoo 
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NcB 



''' -"^max ^ 



(1.3 Mpc) 



(38) 



NcB 



6x1050 



10- 



,SN 



1 cm -^ 



70 1/(3V3) 

103 Arans 



where we have distinguished the average density n^ , close 
to the parent SN (that determines Rmax) from the density 
np in the superbubble or the outer galaxy. The function of 
interest, 7(t), is the real root of the above cubic equation, 
that is: 



7 = 7(7o,^,Xoo;t) 



C 



B 



6*2 + c ei"^ 



c 



B2 



BVB2+4( 
6ct 



1/3 



1 302 

B = -^ + + T. r- 

7o 70 (l + z)x, 



(39) 



The distance travelled by the CB is given by directly 
integrating Eq. (p4|): 



x(7) = Xo 



1 

70 



(40) 



The characteristic distance over which the Lorentz factor 
evolves from 70 to 7o/2 is Xoo/70: roughly 1.3 kpc for 
7o — 10^ and the reference value of Xqo in Eq. (p8|). 

In Fig. (h) we show, for 70 = 10^ and various viewing 
angles 9, the AG flux predicted by Eq. (31), with n = 4.1 
and 7(t) as in Eq. ( |39| ) with the reference value of Xqo 
in Eq. (|38|). For 9^0, and particularly for sufficiently 
large 7 9o, the AG curve described by Eqs. (|3|, |3|) shows 
a very interesting behaviour. Since 7(t) is a decreasing 
function of time, the Doppler factor first increases with 
time, reaches a maximum value at 7(1) ~ 1 and then 
declines. An observer initially outside the beaming cone 
(70 9 > 1), sees an AG that initially rises with time. As 7 
decreases the cone broadens, and around j9^1 beaming 
becomes less efficient, the AG declines. 

In Fig. (^ we show the AG flux predicted by Eq. (^, 
with n = 4.1 and 7(1) as in Eqs. ( pSJ , p^ ), for 70 = 1/9 — 
10^, Up = 10^3 cm^3^ and for various values of the den- 
sity Up^ close to the progenitor SN. For the smaller n?^ , 
the hmiting CB's radius Rmax of Eq. (|l6|) is larger. Conse- 
quently the CB, subsequently ploughing through the ISM, 
loses momentum at a faster pace. The figure shows that 
this may extinguish the AG very soon after the GRB, 
which would make it much harder to observe. This effect 
can also be produced by an increase in the ISM density 
np, relative to the reference choice in Eq. (p^), so that the 
AGs of events not occurring in low-density superbubbles 
would also be hard to observe. These may be (along with 
extinction) the reasons why, in some 50% of cases, GRBs 
appear not to have afterglows. 

For late times, when 76* ^ 1, Eq. ( pq ) implies that 
7 oc t~^/3_ According to Eq. (|3l|), then, the AG light curve 



approaches: 

F^(t) ex t-^; 

r = 1-Fa~ 2.1, 



(41) 



while the i^ dependence stays put at 1/ " . This is compati- 
ble with what is seen in various late-time AG observations 
(see e.g. GRB 980326: Bloom et al. 1999a; GRB 980519: 
Halpern et al. 1999; GRB 990123: Holland et al. 2000a; 
GRB 990510: Stanek et al. 1999, Harrison et al. 1999, Hol- 
land et al. 2000a; GRB 991208: Castro-Tirado et al. 2001); 
GRB 000301c: Masetti et al. 2000, Jensen et al. 2000; GRB 
000926: Fynbo et al. 2001, Harrison et al. 2001, Price 
et al. 2001; GRB 010222: Masetti et al 2001, Stanek et 
al. 2001, Cowsik et al. 2001a). 

11. Comparison with optical observations 

We compare our predictions with raw AG observations, i.e. 
observations not corrected for extinction. The first step in 
our procedure is to work out what a raw ansatz-standard- 
candle supernova, SN1998bw, would look like at the loca- 
tion of each GRB. For that, we use the bare (unextinct) 
SN1998bw deduced by Galama et al. (1998a), we trans- 
port it to the GRB location by way of Eq. (^ and correct 
it for extinction in the host galaxy and in ours. For the 
extinction in our galaxy we use the estimates of Schlegel 
et al. (1998). For the correction at the host (at the emitted 
frequency (1 + z) i') we use the wavelength-dependent ex- 
tinction estimated by the observers from the AG spectrum 
and the colours of the host galaxy. The total extinction of 
the SN contribution for each GRB's optical AG is given 
in Table II. 

The next step in our procedure is to fit the raw data 
minus the SN contribution to a fixed host galaxy luminos- 
ity plus the CB's afterglow. We do not correct the latter 
for extinction, which only affects its fitted normalization. 
In the fits, the afterglow's spectral energy density is given 
by Eq. (31) with R = R,nax and a constant ISM density, 
that is: 



F.-F[7(t)]^"[5(t)] 



l3+Q 



(42) 



with a normalization factor F which is one of the fitted 
parameters, and with 7(1) as in Eqs. (|3^)q The contribu- 
tion of the host galaxy is fixed inside the ±1<t error range 
of the photometry measurements at the late times when 
the CB and the SN have become sufficiently dim. 

The parameters to be fit are F and a in Eq. (B3), as 
well as Xoo, 7o and 9 entering the expressions for d and 
7 as functions of time, Eq. (p9). The fit is done with the 
program MINUIT, checked over and over (in the hunt of 
false local minima) with different input parameters. The 
values of 70, a, 6* and Xoo for the different GRB afterglows 
are listed in Table 11, the overall normalization F will be 



® Admittedly, our expression for the AG light curve is not as 
simple as Eq. (|3|), but it is also analytical. And it is justified. 
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discussed separately. The results of these fits, which are 
very good, are shown in Figs. (R) to (p3). Most of them re- 
fer to R-band observations, which arc the most extensive 
and accurate ones in the optical band, and extend to very 
late times. In order to test that our best fitted parame- 
ters are independent of frequency we have also fitted other 
bands, when sufficiently accurate data are available, and 
obtained very similar best fitted parameters. This can be 
seen in Table II for GRBs 990510 and 990712, for which we 
also present V-band fits. Before discussing the results in 
much more detail in section 13, we pay attention to two 
very peculiar afterglows, and one particularly well mea- 
sured one. 



11.1. GRB 970508, CBs exiting a superbubble? 

The optical AG of GRB 970508 is the only one so far that 
has been seen to rise and fall very significantly. In Fig. ( P3[) 
we show how miserably a fit to this GRB fails, if it is made 
in the same way as all of our other fits. 

We have argued in Section 10.1 that GRB progenitors 
are presumably located in super-bubbles of 0.1 to 0.5 kpc 
size. There may be instances in which the jet of CBs, after 
travelling for such a distance, does not continue onwards 
to a similarly low-density halo region, but encounters a 
higher-density domain. To test whether this may explain 
the very peculiar shape of this AG, we have made a fit with 
two values of the ISM particle-number density, instead of 
one, and a time (or distance from the progenitor) at which 
the transition occurs. The result is shown in Fig. (e3) and 
it is fairly satisfactory. 

The fit parameters correspond to a density increas- 
ing by a factor of ~ 2.2 at t ~ 1.1 day after burst, at 
which point the CBs have travelled some ~ 0.24 kpc, a 
very reasonable radius for a superbubble. The remaining 
parameters are in the usual range, but for 0, which, at 
~ 3.5 mrad, is on the large side. Given this large value 
and 7o ^ 1123, the time at which 7(t)0 = 1 is reached 
exceptionally late; this explains the rise and fall of the 
theoretical curve; see Dar and Dc Riijula (2000a) for an 
earlier version of this result. The relatively large 6 is also 
in accordance with the fact that, in spite of a relatively 
large 7, the equivalent spherical energy of this GRB is 
particularly low, see Tables I and II. 

In an earlier version of this paper, we attributed the 
shape of the AG of GRB 970508 to the effects of gravi- 
tational lensing by an intervening star or binary, of mass 
~ 2 Mq. That was an error. The required mass for an ob- 
ject placed mid-way to the GRB location is almost three 
orders of magnitude bigger. An effect of the observed size 
and shape could also be due to an even heavier object 
(such as a globular cluster) kiloparsecs away from the 
source, but the chance probability for that is negligible. 
The possibility of lensing by stars — which has a few per- 
cent probability and would produce amplification effects 



typically lasting ^^ 1 hour — is still interesting. We discuss 
it in Appendix III. 

11.2. GRB 980425, a very special case? 

As reported in Table I, this GRB is by far the closest and 
yet, its measured 7-ray fiuence is not large. We have ar- 
gued (DD2000a) that in the CB model this is simply due 
to the fact that it is observed at a very large 9: the GRB 
fiuence has the same angular dependence as Eq. (Q). The 
association of this GRB with SN1998bw is clear. In fact, 
the "afterglow" is dominated by the SN. Only the very last 
measured AG point is significantly above the ^^Co decay 
trend of the SN ejecta and would be due to the proper af- 
terglow: that of the CBs (see Fig. 4 of D2000a). With only 
one point above the SN "background" we cannot make, in 
the case of GRB 980425, a detailed fit of the sort we have 
made for the other GRB afterglows. But, as we shall see 
after we gain confidence on the success of the CB model 
in describing X-ray afterglows, we can use the X-ray mea- 
surements for this GRB to determine its parameters and 
to show that, after all, it is not at all a very special case. 

11.3. GRB 990123, the early optical data 

In the case of this GRB, there are good optical data start- 
ing exceptionally early: during the 7-ray activity at t = 
22.18 s after its detected beginning (Akerlof et al. 1999). 
The AG rises abruptly to a second point at t = 47.38 
s, and decreases thereafter. At the earlier stage, the CB 
is still hot and fully ionized, its thermal bremsstrahlung 
(free-free) self absorption is very large and fast decreas- 
ing, resulting in a fast rising AG that turns into a de- 
clining light curve when the CBs become transparent to 
short radio waves (corresponding to optical light in the 
observer frame) . It is possible to explain the initial rise in 
detail, but the scarcity of the data and the surplus of pa- 
rameters make the exercise moot. We choose to describe 
this AG from t = 47.38 s - 5.5 x 10""^ d onwards, the 
first point shown of the measured decline and of Fig. (p5[). 
We can use the relation between local distance and ob- 



server's time, Eq. ( |35[ ), and the specific values of z, 70 and 
60 reported in the Tables for this GRB, to conclude that 
at the start of the optical AG data the CBs are a mere 
0.46 pc away from the progenitor star. This is precisely in 
the domain where the density profile ought to be that of 
Eq. (|34|), n oc r~^, induced by the parent-star's wind and 
ejecta. Since at these early times the CB's deceleration is 
negligible, an r~^ density profile translates directly into 
an optical AG that declines as t~^, see Eq. (pi]). 

In Fig. (E3) we report the result of a fit to the AG 
that includes a term proportional to t~^. The parameters 
of the late (t > 1 d) AG remain essentially unchanged 
relative to the ones used before in constructing Fig. (pi]). 
The normalization of the fitted t~^ term is very close to 
that implied by Eqs. ( pl| ) and (p4|): the first point of the 
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data is exactly reproduced for a density n = 0.54 cm~'^ at 
that point, at which x ~ 0.46 pc. Thus, the CB model also 
succeeds in describing this very early AG in magnitude 
and shapeQ 



12. Comparison with X-ray observations 

The CBs enter their AG phase when they become trans- 
parent to radiation. As we have seen, their X-ray AG is 
dominated by thermal bremsstrahlung and first declines 
with time as ~ 



as described by Eqs. (24, 25). In a 



matter of observer's minutes, for typical parameters, the 
CB's radius reaches its limiting value and the synchrotron 
radiation of Eq. (^, which is proportional to Up R^, be- 
comes the dominant AG source in both the X-ray and op- 
tical domains: the corresponding lightcurves are approx- 
imately proportional. Once again for typical parameters, 
a CB reaches the galactic halo in just a few hours. The 
light curves before that time would be hard to model in 
minute detail, since the ISM density is no doubt changing 
rapidly. Coincidentally, there are not enough continuous 
X-ray measurements during the first few hours after the 
GRBs to reliably extract a density profile. 

The previous considerations justify a very simple de- 
scription of the X-ray light curves. Let W be the (con- 
stant) ratio of X-ray to R-band optical AGs. We expect: 



Fx(t) — Fx(ttrans) 



WFR(t), 



(43) 



where t is the observer's time since the ejection of the 
(last) CB, and FR(t) is the spectral energy density in the 
R-band, i.e., Eq. (|3l|), that we have fitted to the R-band 
AG observations. The two parameters we fit to X-ray light 
curves are Fx(ttrans) x [ttrans]^ and W, resulting in a very 
good description of the measured X-ray afterglows of the 
GRBs of known redshift: 970228, 970828, 971214, 980613, 
990123, 990510, 000926 and 010222. The cases of GRBs 
970508 and 991216 require a non-constant density profile. 
The fits to the X-ray AGs are shown in Figs. (||Q. GRB 
980425 is discussed separately below. 

During a GRB the X-ray luminosity fluctuates as the 
7-ray luminosity does, changing abruptly at the end of 
the GRB into a very fast decline. This is expected in the 
CB model, in which the two behaviours have slightly dif- 
ferent origins: thermal bremsstrahlung from the various 
CBs' surfaces during the GRB, thermal bremsstrahlung 
from the rapidly cooling CBs' volumes as they become 
transparent to the radiation they enclose, in a short time 
ttrans — of the Order of seconds — at the end of each indi- 
vidual CB. In Figs. ( p6| - |3^ ) we have therefore shown our 
fits to data beginning at the start of the sharp X-ray de- 
cline. 



"^ Even the most adamant defenders of fireballs admit that, 
in their scenarios, the absence of "windy" AGs is a problem, 
see, e.g. Piran (2001). 



The use in Eq. (ES) of a more careful treatment of the 
bremsstrahlung contribution (evolving from t~^ to t~^) is 
unwarranted: the assumption of a constant ISM density 
should be inappropriate between ~ 2 x 10~^ and ~ 0.2 
days, and there are no data in that domain except for 
GRB 991216 and perhaps 970508 that, like the early op- 
tical AG of GRB 990123, suggest an initial density vari- 
ation ex 1/r^, resulting in an observed ^ t~^ decline. In 
these cases, we have fitted the data with Eq. (E^) for a 
~ 1/r^ plus constant density along the CB trajectory. For 
the particularly interesting case of GRB 970508, shown 
in Fig. (|^, we have also subtracted from the data the 
contribution of the X-ray line observed before 0.8 d. This 
reduces the two points observed between 0.2 and 1 day by 
a factor ^ 0.39. The overall result is compatible with an 
effect that, at late times, is achromatic, since both the late 
optical and X-ray AGs are proportional to ne; see Fig. (24 ) 
for the optical counterpart. 

In what follows we refer to the initially rapidly-falling 
part of an X-ray AG, as "early" and to the subsequent 
much flatter AG as "late" . 

We have chosen to write the prediction for X-rays as 
in Eq. ( [4^ ) to better expose the expected achromaticity 
of the late AG. Moreover, the data for late X-ray AGs is 
generally more sparse than for optical light; it is therefore 
advisable to exploit the fact that the CB's parameters are 
better fit from the R-band data, and to write the late 
X-ray AG as a rescaling of that data. But the effects of 
absorption are much less severe for X rays and there is 
one parameter, the overall normalization, for which it is 
preferable to use the X-ray fluence in testing the model. 
The values of a, 9, 70, and Xoo, fit for each GRB to the 
R-band AG, are sufficient to deduce, via Eq. (|3l|), the 
shape — but not the normalization — of the expected late 
X-ray fiux in the 2-10 keV domain. This fiux depends, 
via Ue Rmax ^t fixed Xoo , on the number of CBs and their 
individual baryon number, as F ex n^^^ Ncb = Njet (the 
dependence on Bj_, via Vq~^, is extremely weak). Let q, 
for a given GRB, be the ratio of the observed flux to the 
one expected for our canonical Ncb = 6 x 10^°, and for one 
dominant (largest 7) cannonball: n^g — 1. The values of q 
are reported in Table III. They range from ^ 1/3 to ^ 3, 
indicating that the CB model satisfactorily explains the 
late X-ray AG normalizations and that the total baryon 
number of the ensemble of CBs that dominate the AG 
appears to be quite constant. 

One can see in Figs. (|2^) to (|3^), that the X-ray fiu- 
ences at the start of the X-ray decline are ^ lO^'* — 
10~^ ergcm"^ s~^. These figures for the beginning of the 
early X-ray AG compare quite favourably with the typical 
prediction quoted after Eq. (p5|): 4.35 x 10^^ ergcm^^ s^^ 
for a single CB. Once again, the parameters extracted 
from fitting the optical AGs are not sufficient to fix case 
by case the overall early X-ray normalization in Eq. (Ea) , 
which also depends directly on other parameters (notably 
Ncb, ncB and the radii of the still-growing CBs). Using 
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these degrees of freedom we could fit not only the late, 
but also the early absolute X-ray fluences. Suffice it to say 
that the magnitude of the early X-ray afterglow is also the 
one expected in the CB model. 

The fitted values of W in Eq. (||) are in 
fair agreement with the absorption-dependent expecta- 
tion [A(i/x)/A(z/o)] {i^x/i^o)~"°^ from the spectrum of 
Eq. (pO|). Some examples of values of a^^ extracted this 
way are 1.06 ± 0.12 for GRB 970228 (Frontera et al., 
1998), 1.12 ±0.07 for GRB 970508 (Galama et al. 1998b), 
0.95±0.1 for GRB 971214 (Dal Fiume et al. 2000), > 1 for 
GRB 980703 (Vreeswijk et al. 1999) 0.96 ± 0.26 for GRB 
990510 (Pian et al. 2001), 0.9 to 1.1 for GRB 000926 (Piro 
et al. 2001) and 0.97 ±0.05 for GRB 010222 (In 't Zand et 
al. 2001). However, the inferred values of aox are affected 
by large and very uncertain extinctions in the CB and host 
galaxies. Evidence for large extinction of optical AGs of 
GRBs in their host galaxies is provided by the large col- 
umn densities (Nh > 10^^ cm~^) extracted from the X-ray 
observations of some GRBs (e.g., GRBs 970228, 970508, 
970828, 971214, 980329, 980519: Owens et al. 1998; GRB 
980703: Vreeswijk et al. 1999) and from the absorbed spec- 
tra of the optical AG of other GRBs (e.g. GRB 990712: 
Vreeswijk et al. 2000; GRB 991216: Halpern et al. 2000a; 
GRB 000131: Andersen et al. 2000). In fact, the failure to 
detect the optical AG of many long duration GRBs with 
well locafized X-ray AG —like GRB 970111, GRB 970616, 
GRB 970815, and 970828— may be due to strong extinc- 
tion of their optical AG in the host galaxy (Djorgovski et 
al. 2001 and references therein). 



Table II. The data are not very precise and the fit (for 
which we assumed a — 1.1, as fitted for all other GRBs) 
is not one of our best, but it inescapably requires a very 
large viewing angle 9. The best fitted angle is 6* ~ 8.3 
mrad, corresponding, for the fitted 7 ~ 750, to an initial 
Doppler factor S - 37. If the CBs of GRB 980425 had 
been viewed from a typical viewing angle, 9 < I/70, the 
equivalent isotropic energy would have been in the range 
7.3 X 10^^ to 5.8 X 10^2 erg, like that of all other GRBs. 

If we assume that for GRB 980425 the extinction, 
ISM density, and CB radius were the same as for other 
GRBs well measured in X-rays, such as GRB 990510 or 
GRB 010222, we can use Eq. (^l]) to derive the expected 
intensity of the X-ray AG plateau of GRB 980425, see 
Fig. (HI), and its caption. The results are Fx[425] — 
0.32Fx[510], and Fx[425] = 0.15Fx[222], both yielding 
Fx[425] '^ 4 X 10~^^ergcm~^ s~^, in agreement with the 
BcppoSAX observations (Pian, 1999; Pian et al. 2000). 
The double success in deducing a "normal" GRB equiva- 
lent isotropic energy, and the intensity of the X-ray AG, 
constitutes a very strong support for the alleged associa- 
tion of SN1998bw with (a not exceptional) GRB 980425. 

The fitted parameters of the X-ray AG can be used 
to predict the magnitude and shape of the optical AG of 
the blended SN 1998bw/GRB 980425 system, if we assume 
the same V/X extinction ratio as in GRBs 990510, 000926 
and 010222. This we do in the lower part of Fig. (|6|). 
The CBs' contribution dominates at very late time and, 
remarkably, it is in perfect agreement with the HST ob- 
servation (Fynbo et al. 2000) on day 778 after the GRB. 



12.1. GRB 980425 is not exceptional 

The 7-ray fluence of this GRB is not atypical but its red- 
shift, z = 0.0085, is extremely small. If it is not intrinsi- 
cally exceptional, its CBs must be viewed from an atyp- 
ically large angle (DD2000a). For large 9 the CBs' after- 
glow is strongly reduced, as can be seen from Fig. (^, 
allowing for the possibility that the AG is dominated 
by the SN. Consequently, the CB parameters cannot be 
derived from the optical light curve of the blended SN 
1998bw/GRB 980425 system. But they can be deduced 
from the X-ray emission of the system if we assume, unlike 
the observers do (e.g., Pian et al. 1999; Pian et al. 2001) 
that it was produced by the CBs and not by the conven- 
tional quasi-spherical SN ejecta^. Indeed, significant X-ray 
emission from SNe has been detected only at much later 
times after the event. Moreover, the exceptionally slow 
decline of the X-ray AG in this GRB is what is expected 
from the large viewing angle interpretation, see Fig. (0). 
Given all of the above, for this GRB we have "reversed" 
our procedure by first fitting the X-ray AG of the SN 
1998bw/GRB 980425 system. The fit is shown in the upper 
part of Fig. (pq) and the fitted parameters are listed in 



13. Discussion of the results on optical AGs 

With the exception of the AG of GRB 970508, which has 
the sharp "break up" that we have explained via a sudden 
change in density, a look at Figs. (R) to (24) clearly reveals 
that the observed AGs have absolutely no "breaks" . In the 
CB model, the gradual evolution of the proper afterglow 
(that of the CBs) is simply a consequence of the gradual 
decrease of the Lorentz factor 7(t). We give in Table II 
the list of the parameters resulting from our fits (a, 70, 9 
and Xoo) to optical R-band afterglows (and in two cases, 
to the V-band afterglow as well). 

The fits to the CB model are satisfactory, particularly 
since the best-fit parameters turn out to be precisely in the 
expected ranges. On close inspection one notices that our 
curves occasionally undershoot or overshoot some points 
by a small factor, as in GRBs 990123 and 000301c5 This 
is not a surprise: the AG fluences are proportional to the 
ISM number density, which we do not expect to be exactly 
constant for kpc distances, even in the halo of galaxies. If 
such "defects" were not present in our fits, we would have 



We are indebted to E. Pian for discussions on this point. 



^ The feature at t ~ 4 days in GRB 000301c has been in- 
terpreted by Garnavich et al. (2000b) as due to gravitational 
lensing. 
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concluded that the data had been over-parametrized. For 
the same reasons, and because of the systematic errors in 
the data, the values of the parameters we extract from 
our fits should not be taken entirely at face value, even 
though the minimization procedure — which attributes to 
the errors a counterfactual purely statistical origin — re- 
sults in tiny 1 u spreads for the fitted parameters, and in 
X^ values that are in most cases extremely satisfactory. 

All the figures (Q) to ( p4|) refer to optical data for t > 
0.1 days, for which it is reasonable to approximate the 
ISM density by a constant value, describing the density of 
the superbubble and/or the galactic halo. We have already 
discussed the early observations of GRB 990123, for which 
this approximation breaks down. 



13.1. The distribution of fitted parameters 



In the CB model, the parameter a of Eqs. (pQ, 31) is the 
only one for which we have no reason to expect a range of 
different values. It is therefore extremely satisfactory that 
the fitted values of a are, within errors, compatible with 
all of the GRBs having a universal behaviour with the the- 
oretically predicted value: a « 1.1, Eq. (30). The narrow 
distribution of best fitted a values is shown in Fig. ( p7|) . 
In the CB model, we have extracted the values of a from 
the temporal shape of the AG and — adding consistency 
to the picture — they agree well with the values obtained 
from spectral observations, either in X-rays (with spectra 
modified by a best-fit hydrogen column density) or in the 
optical domain (with a galactic colour-dependent extinc- 
tion). Some examples are: 

GRB 970228: ax = 1.06 ± 0.12 (Costa et al. 1997). 
GRB 970508: ao = 1.12 ± 0.07 (Galama et al. 1998c); 
ax = 1.11 ± 0.06 (Galama et al. 1998c). 
GRB 990123: ao = 1.29 ± 0.23 (Holland et al. 1999a). 
GRB 990510: ax = 0.96 ± 0.26 (Pian et al. 2001); 
ao = 1.26 ± 0.15 (Stanek et al. 1999). 
GRB 991208: ao = 1.05±0.09 (Castro-Tirado et al. 2000). 
GRB 991208: ao ~ 1.1 (Takeshima et al. 1999). 
GRB 000301c: ao = 1.15 ± 0.26 (Jensen et al. 2000). 
GRB 000926: ax = 1.2 ± 0.3 (Piro et al. 2001); 
ao = 1.02 ± 0.02 (Sagar et al. 2001a). 
GRB 010222: ax = 0.97 ± 0.15 (In 't Zand et al. 2001); 
ao = 1.07 ±0.09 (Stanek et al. 2001). 

There are also cases for which the reported value of 
a differs significantly from 1.1. One notable instance is 
GRB 990510, for which Beuerman et al. (1999) report a -- 
0.55 and Stanek et al. (1999) find a = 0.61 ± 0.12. To 
extract this value the authors extrapolate the measured 
extinction: E(B - V) = 0.20 (Schlegel et al. 1998). If this 
measured extinction is used to correct only the measured 
B - V values for GRB 990510: 0.57 ± 0.02 (Beuerman 
et al. 1999), and 0.56 ± 0.03 (Stanek et al. 1999), one 
obtains ao = 1.08 ± 0.12. The uncertainties entailed by 
absorption corrections are the reason why we have chosen 



to de-emphasize results that are sensitive to them, whether 
they do, or do not, agree with our expectations. 

The distribution in initial Lorentz factors, 70, shown 
in Fig. (|8|), agrees snugly with our expectation, 70 ~ 10'^, 
extracted from independent information: the fluence and 
the individual-photon energies of GRBs (DD2000a,b) and 
the energies of X-ray lines in their afterglow (DD2001a). 
Notice how surprisingly narrow this distribution is. 

The distribution of viewing angles 6 is shown in 
Fig. (H). The AG data for GRB 000131 consist in only 
three points, while for GRBs 991208 and 000301c the mea- 
surements start rather late. The sensitivity to 6 in our fit 
to these GRBs is not good. We reflect this fact in Fig. (39) 
by having the corresponding entries unshaded. The distri- 
bution is compatible with the expectation that the limited 
experimental sensitivity to GRBs introduces a sharp cut- 
off as 9 increases; see the steep fluence function, Eq. (^. 

The parameter Xqo of Eq. (pq) is the only one for which 
we expect a rather broad distribution. Indeed, it depends 
on the densities close to the GRB progenitor, which ought 
to be quite variable, and in the region where the CB light 
is emitted; see Eq. (|3q). The values of Xqo reported in 
Table II show a spread of a bit over one order of mag- 
nitude, supporting the expectation. In Fig. (BQ) we show 
the distribution of Logio[xoo(Mpc)], which peaks at the 
reference value of Eq. ( p8| ) and extends to smaller values, 
as expected if the average density close to the progenitor 
is sometimes much smaller than our rather large reference 
value: n^^ = 1 cm^'^; and/or the density of the ISM is big- 
ger than our rather low reference value: np — 10^'^ cm^"^. 

The values of a, 6, 70 and Xqo are not sufficient to pre- 
dict the overall normalization of an AG: F in Eq. (p2[), 
whose approximate value is given by the absolutely nor- 
mahzcd Eq. (|3l|). Indeed, F is proportional to the product 
of the number of CBs and their baryon number. To skirt 
absorption corrections we have discussed in section 12 the 
values of F in connection with late X-ray AGs. There, we 
compared the data and the naive expectation for a sin- 
gle (highest-7) dominant CB (n^^ = 1) and our canonical 
NcB , to work out the ratio q of fitted to predicted values 
of F in the X-ray band. The same exercise can be carried 
along for the R-band AGs, with the result that the optical 
values of q are not within a factor of three, but within an 
order of magnitude of q = 1. It is tempting to conclude 
that this may be due to poorly-understood absorption. It 
is in any case clear that the AG magnitude, in the CB 
model, is not a problemFi. 



^° In the fireball model and its descendants, the efficiency of 
conversion of kinetic energy to photons is claimed to be high, 
and both the GRBs and their AGs are due to the same mecha- 
nism: synchrotron radiation. It is therefore difficult to explain 
why, at the end of the GRB, the radiation rate suddenly drops 
by two orders of magnitude or more, and why there is more 
energy in the GRB than in the afterglow (see, e.g. Burenin et 
al. 1999; Pian et al. 2001, In 't Zand et al. 2001). 
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The GRB consists of the photons emitted by the hot 
CBs as they exit the SN shell. In the rest system of a CB, 
the individual energies, as in Eq. (p^, are a fraction (1 + 
z)/(5o of the observed energies, with Sq = S[^o,6]. The total 
energy emitted by a GRB, in the rest system of its CBs, 
is in the form of isotropically distributed photons that 
appear to us as the coUimated 7-ray burst. This comoving 
total energy is related to the observed fluence F^ by: 



E 



CB 






(44) 



For the GRBs all of whose parameters are well determined, 
we list in Table III the values of E^;^ . Interestingly, their 
distribution — shown in Fig (ETI) — is also quite narrow. 

We have not discussed in this subsection the parame- 
ters of GRB 980425, hsted in Table II. They are obtained 
from a fit to the X-ray — as opposed to optical — AG, and 
they are imprecise. The deduced value of E5^^ is 0.16 x 10"''* 
erg, a bit lower than those listed in Table II. This is to be 
expected, the small S of this GRB makes its GRB softer, 
and less prominent within the BATSE energy window. 

To summarize, the distributions of parameters are in 
extremely good agreement with the expectations of the 
CB model and, if anything, they are astonishingly close 
to what they would be for "standard candle" GRBs. 

13.2. The GRB/SN association in view of our results 

It is useful to discuss the evidence for a SN component in 
the GRB optical AGs in order of decreasing redshift. The 
fact that we have a consistent and successful description 
of optical afterglows strengthens the interpretation of this 
putative evidence. 

Examining Figs. (0) to (|24|), we draw the following 
conclusions. In the six more distant GRBs, ranging from 
GRB 000131 at z = 4.5 to GRB 010222 at z> 1.474, there 
is no evidence for or against a SN 1998bw-like component. 
In GRB 000418, at z = 1.119, there is an indication of 
an excess, compatible with the SN. In GRB 980613, at 
z = 1.096, the evidence, though based mainly on just one 
point, is very strong. In GRB 991216, at z = 1.02, there is 
a clear indication of a late excess over the CB's afterglow, 
though a SN1998bw-hke contribution does not describe it 
very well (a slightly earlier bump would do a very good 
job, indicating that the standard-candle hypothesis for the 
SN is good, but not perfect). In GRB 980703, at z = 0.966, 
the SN excess is visible and well fitted, but the errors are 
large. For GRB 970828, at z = 0.957, there are no AG 
observations at optical wavelengths. In the peculiar case 
of GRB 970508, at z = 0.835, there is in the data a clear 
excess at late times that is very well fitted by our SN 
ansatz. For the next three closer GRBs (991208, 970228, 
990712), at z = 0.706 to 0.434, the evidence is completely 
convincing that a SN1998bw-like contribution is required 
to fit the data. In the case of GRB 990712, once again, a 



SN peak occurring slightly earlier than that of a redshift- 
corrected SN1998bw would provide a better description 
of the AG. Finally, GRB 980425, at z = 0.0085 is indeed 
associated to a SN: our fairly satisfactory standard-candle 
choice. 

A clear trend is apparent in the last paragraph. The 
closer a GRB is, the better the evidence for its association 
with a SN. The trend is entirely consistent with the fact 
that, for the more distant GRBs, a SN contribution to 
the AG could not be seen, even if it was there. In all 
cases where the SN could be seen, it was seen, with the 
evidence gaining in significance as the distance diminishes. 
The temptation to conclude that all long-duration GRBs 
are associated with SNe appears to us to be irresistible, 
even if an irrefutable proof will never be possible. 

13.3. Do cannonballs deserve their name? 

We have argued that CBs should reach an asymptotic ra- 
dius in a very short time and travel thereafter as literal, 
i.e. non-expanding, "cannonballs" . Can this statement be 
contrasted with the data? To answer this question we have 
analyzed the AGs produced by CBs whose radius inertially 
increases at a fixed speed, /3cxp c, in their rest system. The 
details are given in Appendix II. The result is that the late 
AGs behave in this case as F^ oc t"'^, with r = 9 {l + a)/5. 
For a ~ 1.1, as we have argued, r ~ 3.8, which completely 
disagrees with the AG light curves. For t == 2.1, in agree- 
ment with the latter, a ~ 0.17, which completely disagrees 
with the measured spectra. Thus, the cannonballs of the 
CB model do deserve their name. 



14. Summary and conclusions 

We have previously argued that the cannonball model of- 
fers a successful and simple explanation of the fiuence, 
energy spectrum, and temporal behaviour of the prompt 
7-rays of a GRB (DD2000b). ^From these considerations 
we extracted a CB's typical Lorentz factor, 70 ^ 10'^, and 
typical baryon number Ncb^ in the vicinity of 6 x 10^°, or 
0(10) times as much for a jet of CBs in a multipulse GRB 
(DD2000b, DD2001b). Using only these CB-related input 
parameters and a reasonable initial CB expansion veloc- 
ity, we have explicitly worked out all of the properties of 
X-ray and optical afterglows. 

As an intermediate result, we derived the temporal be- 
haviour of the radius of a CB and showed that, faithful 
to their name, CB radii tend to a constant Rmax in mere 
minutes of observer's time. The value of Rmax depends on 
the ISM density close to the GRB progenitor. GRBs whose 
CBs have large radii have afterglows whose temporal de- 
cline is very fast. This is also the case for GRBs whose CBs 
travel through a relatively dense ISM. These may be the 
reasons, along with strong extinction in the host galaxy, 
why the search for AGs is not always successful. 
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We have shown how well the CB model describes 
all the properties of the X-ray and optical afterglows of 
GRBs. Our results do not — and could not — take into ac- 
count possible variations of the ISM density along a CB's 
path; they are in this sense "descriptions" rather than fits. 
In spite of this, the descriptions are excellent and the con- 
sistency of the results is impressive. 

The observed behaviour of both X-ray and optical AGs 
is the predicted one. All the parameters extracted from the 
fits have values or distributions close to the expected ones. 
This is the case for the integrated fluences in the early and 
late X-ray and optical bands, for the values of the Lorentz 
factor 7, for the distribution of observation angles 0, for 
the spectral index a (that we determine, not from the 
spectra themselves, but from the time-dependence of the 
late afterglow) , and for the parameter Xqo of Eq. ( |3q ) , that 
governs the pace of the CB's slowdown, and is a combi- 
nation of the ISM densities, and the CB's baryon number 
and radius. Throughout, we have chosen to de-emphasize 
the results that are most sensitive to systematic errors, 
such as absorption corrections. Thus, we have not system- 
atically extracted parameters from flux normalizations, 
nor reported any x^ tests of the quality of the fits (which 
would be misleadingly good) or statistical-error estimates 
on the fitted parameters (which would be misleadingly 
small) . 

The distribution of the derived quantity E5;^ (the en- 
ergy of the GRB photons in the CBs' rest system) has a 
range of less than a factor of four, making long-duration 
GRBs not far from standard candles. We have demon- 
strated that GRB afterglows can be understood in detail. 
Perhaps this will pave the way to the use of GRBs as 
cosmological beacons. 

GRB 980425 is at first sight a special case: its 7-rays 
are rather soft, it is extremely near by cosmological stan- 
dards, and very clearly associated with a SN. In the CB 
model it is not exceptional, only seen at a relatively large 
angle. Its X-ray afterglow is normal and emitted by its 
CBs, not by the isotropic ejecta of SN1998bw. Its optical 
AG is dominated by the SN for almost two years, but its 
last measured point is due to the CBs, and agrees with 
the expectation. 

Very early optical AG data are only available for GRB 
990123. The CB model is also capable of describing them 
naturally: their magnitude and time-dependence are those 
expected for CBs moving through the density profile pro- 
duced by winds in the Wolf-Rayet phase of the progenitor. 
Future early observations should test this feature of the 
CB model: a linear dependence of the AG fiuence on the 
varying local-density profile. 

For GRB 970508 we can provide a good fit to its pe- 
culiar AG only if we assume that its CBs encountered a 
sudden change in ISM density, as they would if they are 
exiting a superbubble into a higher-density region. 

The very early X-ray and optical AGs are not achro- 
matic: they are dominated by two mechanisms with dif- 



ferent time-dependence: thermal bremsstrahlung and syn- 
chrotron radiation. But in a matter of a few hours all of 
the AG is dominated by synchrotron radiation, and the 
light curves should be achromatic, as observed. 

Our descriptions of the optical AGs indicate or require 
the contribution to the light curve of a supernova akin to 
SN1998bw, in all eight out of sixteen cases for which the 
event occurred close enough for such a contribution to be 
observable. The conclusion that there is a roughly one- 
to-one association between core-collapse SNe (or perhaps 
just type lb and Ic SNe) and long-duration GRBs is very 
tempting. The enormous beaming factor of the CB model 
makes this conclusion tenable and consistent. 

We are currently completing the study of the predic- 
tions of the CB model for radio afterglows, for which self- 
absorption in the CBs is relevant, and requires a careful 
analysis. The results, soon to be announced (Dado et al. 
2002), are excellent. We also plan to discuss elsewhere the 
interesting implications of the CB model for cosmic-ray 
physics, and for other accreting compact objects that eject 
relativistic jets, such as radiogalaxies, AGNs, quasars, mi- 
croquasars, blazars and microblazars. 

It would be interesting to compare our results to those 
of conically ejected shocked fireballs. We have argued that, 
for this, it would be more convincing not to have the ob- 
served firetrumpets pointing precisely to planet Earth. For 
the time being, we have proved that the CB model — 
based on simple and definite hypothesis — makes predic- 
tions that are univocal (as opposed to multiple-choice), 
explicit, analytical in fair approximations, quite simple, 
very complete, and very successful. 
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Appendix I: Subdominant AG mechanisms 

Three mechanisms declining less fast than thermal 
bremsstrahlung — but typically subdominant relative to 
synchrotron radiation — contribute to GRB afterglows: 
relativistic bremsstrahlung from Coulomb collisions of CB 
electrons with the ISM constituents, atomic transitions 
of CB atoms excited or ionized by these same collisions, 
and inverse Compton scattering of the CBs' electrons on 
the cosmic background radiation. We discuss them here 
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to substantiate our assertion that thermal bremsstrahlung 
and synchrotron radiation dominate the AGs. 

Relativistic bremsstrahlung 

The electrons of the CB, in their highly relativis- 
tic collisions with the ISM nuclei, emit non-thermal 
bremsstrahlung radiation. The total power per unit area 
received by the observer from this source is: 

dFRB 45.7argNcBnpmec'^7^J^ 

dtdr*"* AttBI ' ^ ^ 

where ro ~ 2.18 fm is the classical radius of the elec- 
tron and a w 1/137 is the fine structure constant. For 
NcB = 6 X 10^", 7 = 5= 10^, z = 1 and Up = Icm^^^ 
the above expression yields 3.8 x 10"^^ erg/(cm^ s). This 
radiation has the same dependence on 7 and 6 as the 
synchrotron radiation of Eq. (p7|), but it is some two 
orders of magnitude smaller. This spectrum extends to 
~me5/(l-|-z) r^ 250 MeV; it may be a contribution to the 
very high energy 7-rays observed by EGRET and COMP- 
TEL in some GRBs. 

Atomic transitions 

The recombination and de-excitation of the CB atoms ex- 
cited and dissociated by Coulomb collisions with the ISM 
particles contributes to the AG with a power: 



1 TT r: Un m, 



c3NcBln[A]7(54 



dFc 
dtdrj 

where In [A] is the Coulomb logarithm, 



inDl 



(46) 



In [A] = In 



^2mec2 7'T„ 
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I = 13.6 eV is the ionization energy of hydrogen, and Tmax 
is the maximum kinetic energy that can be imparted to a 
stationary electron in a single collision with a relativistic 
particle of mass M: 



T — 



2 me C^ 7^ 



l-|-2 72mo/M-|-(mc/M)2 ■ ^ ' 

For NcB = 6 X 10^°, j = 5= 10^, np = 1 cm'^ and z = 1, 
Eq. (E^) yields 7.15x 10""'^^ erg/(cm2 s), which is negligible 
with respect to the synchrotron radiation of Eq. ( p7| ) . This 
power is radiated mostly as synchrotron emission from the 
knocked-on electrons and by line emission from hydrogen 
recombination (boosted to X-ray energies). 

Inverse Compton scattering 

The scattering of the CBs' electrons on ambient starlight 
and the cosmic background radiation (CBR) produces a 
radiated power: 

dFcBR 32^r2p^cNcB7''5^ 



dtdf7 
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(49) 



where p^ is the energy density of the radiation field [p^ — 
0.24 (1 + z)4 eVcm-3 for the CBR). For Ncb = 6 x 10^^ 
7 = (5 = 10'^ and z = 1, the above expression yields 
1.5 X 10~^'* erg/(cm2 s) for the CBR contribution, which 
is negligible relative to synchrotron radiation up to a very 
late afterglow phase, when ne may be very small. The con- 
tribution of Compton scattering to a polarization of the 
signals, on the other hand, may not be negligible (Shaviv 
and Dar, 1995). 

Inverse Compton scattering of the CBs' high en- 
ergy electrons, with a proper spectrum dug/dE ~ 
EjT^'^, on their self-produced synchrotron and thermal 
bremsstrahlung radiation produces very high energy pho- 
tons with a spectrum dn^/dE ~ E~2i. Por distant GRBs, 
this spectrum is cutoff at sub TeV energies by pair produc- 
tion on the infrared background radiation. But, for very 
nearby GRBs, it may be observable up to extremely high 
energies, larger than those observed from blazars. 

Appendix II: Ever expanding "cannonballs" 

It is very instructive to study the possibility that CBs, 
instead of reaching an asymptotic radius, would continue 
to expand significantly. To find 7(1) in this case, Eq. (36) 
is to be substituted by: 

dx 
dR = /?exp— , (50) 

7 

while Eqs. (p4l p5| ) remain unchanged. Insert Eq. (pw into 
Eq. (M) and integrate, to obtain, for constant Up: 



ix ['J } -^trails ' ^inax' 



1 



1 



(47) R: 



3 _ 



3 Ncb A 



cxp 



2 7rnr, 



(51) 



The CBs reach an asymptotic radius as 7 —f 1, which 
bears some resemblance to that of Eq. (|iq), but it is 
reached much later and it is much larger (e.g. R ~ 
3.8 X 10^^ cm for np = 10"^ cm-^, /3trans = l/(3%/3), 
and Ncb = 6 x 10^°). 

Upon insertion of Eq. 
analogue of Eq. (^ : 



l|) into Eq. (^4]), we obtain the 



70 



d7 



1 



92^2 



ry8/3 

Ncb 






-2/3 



2ct 



7rR2 



3(1 +z)Xoo 

(52) 

The integral can be 
done analytically, but is not compact. 

For late times, when 76* ^ 1, Eq. ( p2| ) implies that 
7 oc t~'^/^. According to Eq. (pl|), then, the AG light curve 
approaches: 



where we have neglected R^j.. 



JK 



F^(t) ex t- 



^ (l + a)~3.8 
5 



(53) 
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This behaviour cannot be reconciled with the data, as ex- 
plained in the text. 

Appendix III: gravitational lensing of moving CBs 

The phenomenon of gravitational lensing is well known. A 
lensing object of mass M has a Schwarzschild radius Rs = 
2GnM/c^. If Da(z) is the angular distance from source 
to observer and x is the fractional distance to the lens, the 
Einstein radius of the system is Re = [2Rs Dax(1— x)]^/^. 
As the lensing object crosses close to the line of sight (or, 
in our case, as the line of sight to the fast-moving CB 
passes close to the lensing object) the amplification A is: 



A rough estimate of e taking into account that stars 
gather in galaxies gives a similar result. Let the surface- 
number density of stars in a galaxy, as a function of 
distance to the centre, be approximated by S*(r) = 
S*(0) e^''/'^, with h^ 5 kpc. For a reference galaxy with 
N, = 10" stars i;(0) ~ 640 pc^^. Define a galaxy's ef- 
fective lensing radius so that S»(reff)x^ Re ~ 1. For the 
quoted values of xj_ and Re this means reff ~ 9 kpc. 
Approximate the surface density of galaxies at z < 1 by 
the observed value for galaxies with R-magnitude below 
25: Eg ~ 4.6 x 10^ rad'^ (Casertano et al. 2000). The 
lensing probability at an angular distance Da is then 



e^7rr2g.EG/Di^4%. 





1 Vu + u^ ' 


u(t) = 


^^2 1 1- ^inaxj 


r(t) EE - 


Re 


i(t) ~ c 


1 1 ^ ' 



1/2 



(54) 



where Umin is the minimum distance to the lens, in Ein- 
stein radii, of the line of sight to the CB during its motion. 
Gravitational lensing of a moving CB is peculiar in two 
ways: the apparent velocity is superluminal, and the time 
"width" of the effect, r(t), is itself time dependent, since 
the CB is decelerating as the lensing occurs: 



^(t)=r(0)^=r(0) 

V_L(t) 



70 



1 2 



7(t) 



I + [010? 



(55) 



For a solar-mass star placed halfway to a GRB at 
z = 1, the typical duration of a lensing event is r(0) ^ 1 
hour. The average Einstein radius of a solar-mass star 
placed somewhere on the way to such a location is Re [©] = 
Re (Mq/M)1/2 2 ([x(l - x)]i/2) ^ I860 AU. What is the 
optical depth (or apriori probability), e, for an observ- 
able lensing by such an object? Consider lensing dur- 
ing the first 10 days of an AG, when its fluence is rel- 
atively high and during which, for typical parameters, 
the CBs travelled a (local) distance x ~ 2 kpc. The ap- 
parent transverse distance is x_l = x 6*7(1 -j- z) ~ 1 pc. 
The average luminosity density of the local Universe is 
p ^ (1.8 ± 0.2) h X lO^Lo/Mpc^ and the mass to lumi- 
nosity ratio of star populations is M/L ^ 5 to 10 in solar 
units, so that the number density of "typical" solar-mass 
stars is uq - pM/L - 8.8 x lOVMpc^ for h - 0.65. The 
optical depth is: 



XiR0DAn0((l+z)3) 



(56) 



In the interval extending to z = 1, and for our adopted 
cosmological parameters, the volume average {{1 + z)^) is 
~ 5. Thus we obtain e ^ (4±2)%, which makes the lensing 
effects hopefully visible. 
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Table I - Gamma-ray bursts of known redshift 



Table II - The CB, host-galaxy, and extinction pa- 
rameters 



GRB 


z 


Dl 


F^ 


E^ 


RHG 


970228' 


0.695 


4.55 


1.1 


0.22 


25.2'** 


970508^ 


0.835 


5.70 


0.49 


0.10 


25.0'« 


970828^ 


0.957 


6.74 


9.6 


2.06 


24.52" 


971214"* 


3.418 


32.0 


0.94 


2.11 


25.6^' 


980425^ 


.0085 


.039 


0.44 


8.1E-6 


14.3^2 


980613*^ 


1.096 


7.98 


0.17 


0.61 


24023 


980703^ 


0.966 


6.82 


2.26 


1.05 


22.6^'' 


990123** 


1.600 


12.7 


26.8 


19.80 


23.9^^ 


990510^ 


1.619 


12.9 


6.55 


5.00 


27Q26 


990712'" 


0.434 


2.55 


6.5 


0.53 


21.8^^ 


991208" 


0.70 


4.64 


10.0 


1.51 


24428 


991216'^ 


1.020 


7.30 


19.4 


5.35 


24.82^ 


000131" 


4.500 


44.4 


4.2 


11.60 


27830 


000301c" 


2.040 


17.2 


0.41 


0.46 


28.0^' 


000418'^ 


1.119 


8.18 


2.0 


0.82 


23.9^2 


000926'** 


2.066 


17.4 


2.20 


2.60 


25.6^^ 


010222'^ 


1.474 


11.5 


12.0 


7.80 


25.9^"' 



Comments: z: Redshift. Dl: Luminosity distance in Gpc 
for n^ = 0.3, riA = 0.7 and Ho = eSkms^i Mpc^i. F^ 
BATSE 7-ray fluences in units of 10~^ erg cm~^. E-y 
(Equivalent spherical) energy in units of 10^^ ergs. R[HG] 
R- magnitude of the host galaxy, except for GRB 990510, 
for which the V-magnitude is given, corrected for galactic 
extinction 

References: 1: Djorgovski et al. 1999a; 2: Metzger et 
al. 1997; 3: Djorgovski et al. 2000; 4: Kulkarni et al. 1998b; 
5: Tinney et al. 1998; 6: Djorgovski et al. 1998a; 7: Djor- 
govski et al. 1998b; 8: Kelson et al. 1999; 9: Vreeswijk et 
al. 1999a; 10: Hjorth et al. 1999; 12: Vreeswijk et al. 1999b; 
13: Andersen et al. 2000; 14: Feng et al. 2000; 15: Bloom 
et al. 2001; 16: Fynbo et al. 2001; 17: Jha et al. 2001; 
Fruchter et al. 2001; 18: Fruchter et al. 1999b; 19: Plan 
2001; 20: Djorgovskyet al. 2001 21: Odewahn et al. 1998 
22: Galama et al. 1998a; 23: Djorgovski et al. 2000; 24 
Bloom et al. 1998b; 25: Bloom et al. 1999b; 26: Plan 2001 
27: Hjorth et al. 1999; 28: Diercks et al. 2000; 29: Djor- 
govski et al. 1999b; 30: Andersen et al. 2000; 31: Smette 
et al. 2000; 32: Bloom et al. 2000; 33: Fynbo et al. 2001; 
34: Fruchter et al. 2001. 



GRB 


70 


a 


9 


Xoo 


R[HG] 


AsN 


970228 


540 


1.10 


1.686 


0.155 


25.55 


0.50 


971214 


999 


1.20 


0.708 


0.373 


25.69 


0.94 


980613 


509 


1.09 


1.619 


0.241 


24.07 


0.82 


980703 


779 


1.08 


0.953 


0.344 


22.54 


0.88 


990123 


1325 


1.09 


0.420 


0.954 


23.90 


0.96 


990510 


991 


1.10 


0.261 


0.777 


27.80 


0.50 


990510" 


907 


1.08 


0.318 


0.504 


27.80 


0.40 


990712 


948 


1.09 


0.750 


1.191 


21.93 


0.50 


990712'' 


957 


1.08 


0.863 


1.319 


22.57 


0.40 


991208 


1034 


1.26 


0.100 


1.357 


24.81 


0.80 


991216 


972 


1.09 


0.375 


0.953 


24.64 


0.80 


000131 


1200 


1.26 


0.100 


0.793 


27.80 


0.93 


000301c 


1040 


1.19 


2.223 


0.141 


28.00 


0.90 


000418 


1017 


1.17 


0.970 


1.961 


23.74 


0.92 


000926 


760 


1.20 


0.740 


0.133 


25.63 


0.96 


010222 


1178 


1.10 


0.465 


1.026 


25.76 


0.94 


970508'' 


1123 


1.10 


3.51 


0.293 


24.69 


0.94 


980425*' 


769 


1.10 


8.30 


0.252 


14.30 


0.93 



Comments: 70: Initial Lorentz factor. 9: Viewing angle 
relative to the CB line of motion, in milliradians. Xqo: CB 
slow-down parameter, in Mpc (7 — 7o/2 at x = Xoo/7o)- 
R[HG]: Fitted value of the R-magnitude of the host galaxy 
(except for GRB 990712, for which also the V-magnitude 
is given) not corrected for galactic extinction. Asn: Atten- 
uation of the SN1998bw-like contribution due to galactic 
extinction, a: V-band afterglow parameters, b: Two GRBs 
are special: GRB 970508 is fit with two constant ISM den- 
sities, the Xoo quoted value corresponds to the initial one; 
in GRB 980425 the SN outshines the CBs in the optical, 
the fit is to the X-ray AG, a = 1.1 was assumed, and the 
parameter determinations are very imprecise. 
References: 

GRB 970228: Castander, et al. 1999a; Castander, et 
al. 1999b; Djorgovski et al. 1999a; Fruchter et al. 1997a: 
Galama et al. 1997; Galama et al. 2000; Garcia et al. 1998 
Guarnieri, et al. 1997; Metzger et al. 1997a; Pedichini 
et al. 1997; Sahu et al. 1997a; Sahu et al. 1997b; van 
Paradijs et al. 1998. 

GRB 970508: Bloom et al. 1998b; Castro-Tirado et 
al. 1998b; Chevalier & Ilovaisky 1997; Djorgovski et 
al. 1997; Fruchter et al. 1997b; Fruchter et al. 2000; 
Galama et al. 1998b; Metzger et al. 1997b; Pedersen 
et al. 1998a; Schaefer et al. 1997; Sokolov et al. 1997; 
Sokolov et al. 1998; Zharikov et al. 1998. 
GRB 971214: Diercks et al. 1998; Halpern et al. 1998b; 
Kulkarni et al. 1998b. 

GRB 980613: Djorgovski et al. 1998a; Djorgovski et 
al. 2000; Halpern et al. 1998a; Hjorth et al. 1998. 
GRB 980703: Bloom et al. 1998b; Castro-Tirado et 
al. 1999a; Holland et al. 2000b; Holland et al. 2001; 
Pedersen et al. 1998b; Sokolov et al. 1998; Vreeswijk et 
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al. 1999c; Zapatero Osorio et al. 1998. 

GRB 990123: Castro-Tirado et al. 1999b; Fruchter et 

al. 1999a; Galama et al. 1999a; Holland et al. 2000a; 

Kulkarni et al. 1999. 

GRB 990510: Beuermann et al. 1999; Covino et al. 1999; 

Fruchter et al. 1999c; Galama et al. 1999b; Harrison et 

al. 1999; Holland et al. 2000a; Marconi et al. 1999a,b; 

Pietrzynski & Udalski 1999a,b,c; Stanek et al. 1999. 

GRB 990712: Hjorth et al. 2000a; Sahu, et al. 2000 

GRB 991208: Castro-Tirado et al. 2001; Sagar et 

al. 2000a. 

GRB 991216: Djorgovski et al. 1999b; Garnavich et 

al. 2000a; Halpern et al. 2000a; Sagar et al. 2000a. 

GRB 000131: Andersen et al. 2000. 

GRB 000301c: Jensen et al. 2001; Halpern et al. 2000b; 

Garnavich et al. 2000b; Masetti et al. 2000; Rhoads & 

Fruchter 2001; Sagar et al. 2000b; Veillet 2000a. 

GRB 000418: Berger et al. 2000; Henden et al. 2000; 

Klose et al. 2000; Metzger & Fruchter 2000. 

GRB 000926: Fynbo et al. 2001; Halpern et al. 2000C 

Harrison et al. 2001; Hjorth et al. 2000b; Price et al. 2001; 

Sagar et al. 2001a; Veillet 2000b. 

GRB 010222: Cowsik et al. 2001; Jha et al. 2001; 

Fruchter et al. 2001; Garnavich et al. 2001; Masetti et 

al. 2001; Oksanen et al. 2001; Orosz et al. 2001; Price et 

al. 2001; Sagar et al. 2001b; Stanek et al. 2001; Valentini 

et al. 2001; Watanabe et al. 2001; Veillet 2001. 



Table III - The rest frame GRB energy and X- 
ray AG of GRBs with measured redshift and X- 
ray AG from their observed 7-ray fluence and the 
optical AG parameters 



GRB 


z 


Dl 


F-r 


70 


So 


q 




970228 


0.695 


4.55 


1.10 


540 


591 


0.86 


0.78 


970508 


0.835 


5.70 


1.10 


1123 


137 


1.26 


1.47 


970828 


0.957 


6.74 


1.10 


1153 


1160 


0.77 


1.34 


971214 


3.418 


32.0 


0.94 


999 


1331 


1.28 


1.11 


980613 


1.096 


7.98 


0.17 


509 


606 


1.34 


1.74 


990123 


1.600 


12.7 


26.8 


1325 


2023 


1.45 


1.84 


990510 


1.619 


12.9 


6.55 


991 


1858 


2.60 


0.78 


991216 


1.020 


7.30 


19.4 


972 


1716 


0.38 


1.54 


000926 


2.066 


17.4 


2.20 


761 


1115 


2.12 


0.60 


010222 


1.474 


11.5 


12.0 


1109 


1812 


1.43 


1.31 



Comments: Dl: Luminosity distance, for Q,^ — 
0.3, flA = 0.7 and Hq = 65kms~i Mpc^^, in Gpc. 
F^: BATSE/BeppoSAX 7-ray fluences in units of 10^^ 
erg cm^^. 60: Initial Doppler factor, q: The ratio between 
observed and predicted late-time X-ray fluxes in the 
2-10 keV band, for a single unextinct standard CB with 
NcB = 6 X 10^° . E™: Energy radiated by the ensemble 



of CBs in its rest frame, in units of 10 erg 
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(b) CANNONBALL 

Fig. 1. — (a) A firecone or, more properly, a firetrum- 
pet. In the scenario discussed in the text (initiated by 
Rhoads 1997), the cone expands conically for a distance, 
after which the jet angle 0j widens as its front travels, 
(b) Cannonballs (shown here, somewhat pedantically, a 
bit Lorentz-contracted) subtend decreasing angles as they 
travel. The only relevant angle in the CB model is the 
observer's viewing angle 9. 




Fig. 2. — An "artist's view" (not to scale) of the CB 
model of GRBs and their AGs. A core-collapse SN re- 
sults in a compact object and a fast-rotating torus of non- 
ejected fallen-back material. Matter (not shown) catas- 
trophically accreting into the central object produces a 
narrowly coUimated beam of CBs, of which only some of 
the "northern" ones are depicted. As these CBs pierce the 
SN shell, not precisely on the same spot, they heat and 
re-emit photons, that are Lorentz-boosted and coUimated 
by the CBs' relativistic motion. 
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Fig. 3. — A CB, in its rest system, sees the constituents 
of the ionized ISM impinge in one direction. The CB's 
chaotic magnetic field disperses these particles so that 
they come out isotropically. Electrons, but not protons, 
lose their energy by synchrotron radiation. 
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Fig. 4. — The radius of a CB as a function of observer's 
time, after the CB becomes transparent to radiation, for 
two choices of the initial transverse expansion velocity 

Ptrans C. 
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Fig. 5. — Afterglow flux, in arbitrary units, as a function 
of observer's time, for 70 — 10"^ and various viewing angles 
6, as given by Eqs. ( pi|) with n = 4.1 and 7(t) as in Eq. (39) 
with the reference value of Xoo in Eq. Ifi 
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Fig. 6. — Afterglow flux, in arbitrary units, as a function 
of observer's time, for 70 — 1/0 = 10"^ and various values 
of the (constant) ISM density, n^^, close to the GRB pro- 
genitor, as given by Eqs. (|3l|), with n = 4.1 and 7(t) as in 
Eq. (H). 
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Fig. 7. — Comparisons between our fitted R-band after- 
glow (upper curves) and the observations listed in Table ff , 
not corrected for extinction, for GRB 970228, at z = 0.695. 
Upper panel: Without subtraction of the host galaxy's 
contribution (the straight line). Lower panel: With the 
host galaxy subtracted, and the CB's AG (the line of 
squares) given by Eqs. (|3l|) and (^). The contribution 
from a 1998bw-like supernova placed at the GRB's red- 
shift, Eq. (0), corrected for extinction, is indicated in both 
panels by a line of crosses. The SN bump is clearly dis- 
cernible. 



Fig. 8. — Comparisons between our fitted R-band after- 
glow (upper curves) and the observations listed in Table II, 
not corrected for extinction, for GRB 971214, at z = 3.418. 
Upper panel: Without subtraction of the host galaxy's 
contribution (the straight line). Lower panel: With the 
host galaxy subtracted, and the CB's AG (the line of 
squares) given by Eqs. ( plj ) and (p9|). The contribution 
from a 1998bw-like supernova placed at the GRB's red- 
shift, Eq. (0), corrected for extinction, is indicated by a 
line of crosses. The expected supernova contribution is too 
weak to be observed. 
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Fig. 9. — Comparisons between our fitted R-band after- 
glow (upper curves) and the observations listed in Table II, 
not corrected for extinction, for GRB 980613, at z = 1.096. 
Upper panel: Without subtraction of the host galaxy's 
contribution (the straight line). Lower panel: With the 
host galaxy subtracted, and the CB's AG (the line of 
squares) given by Eqs. (|3l|) and (^). The contribution 
from a 1998bw-like supernova placed at the GRB's red- 
shift, Eq. (0), corrected for extinction, is indicated in both 
panels by a line of crosses. A SN 1998bw-like contribution, 
though based on only one significant point, appears to be 
required. 



Fig. 10. — Comparisons between our fitted R-band 
afterglow (upper curves) and the observations listed in 
Table II, not corrected for extinction, for GRB 980703, 
at z = 0.966. Upper panel: Without subtraction of the 
host galaxy's contribution (the straight line). Lower panel: 
With the host galaxy subtracted, and the CB's AG (the 
line of squares) given by Eqs. ( pi| ) and (|3|). The contri- 
bution from a 1998bw-like supernova placed at the GRB's 
rcdshift, Eq. (|^), corrected for extinction, is indicated by 
a line of crosses. A SN 1998bw-like contribution, though 
the errors are large, appears to be required. 
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Fig. 11. — Comparisons between our fitted R-band after- 
glow (upper curves) and the observations listed in Table II, 
not corrected for extinction, for GRB 990123, at z = 1.6. 
Upper panel: Without subtraction of the host galaxy's 
contribution (the straight line). Lower panel: With the 
host galaxy subtracted, and the CB's AG (the line of 
squares) given by Eqs. (|3l|) and (^). The contribution 
from a 1998bw-like supernova placed at the GRB's red- 
shift, Eq. (0), corrected for extinction, is indicated by a 
line of crosses. The fit is not very good. There are no ob- 
servations at the time when a SN1998bw-like contribution 
may have been seen. 



Fig. 12. — Comparisons between our fitted R-band after- 
glow (upper curves) and the observations listed in Table II, 
not corrected for extinction, for GRB 990510, at z = 1.619. 
Upper panel: Without subtraction of the host galaxy's 
contribution (the straight line). Lower panel: With the 
host galaxy subtracted, and the CB's AG (the line of 
squares) given by Eqs. ( plj ) and (p9|). The contribution 
from a 1998bw-like supernova placed at the GRB's red- 
shift, Eq. (0), corrected for extinction, is indicated in both 
panels by a line of crosses. The expected supernova con- 
tribution is too weak to be observed. 
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Fig. 13. — Comparisons between our fitted V-band after- 
glow (upper curves) and the observations listed in Table II, 
not corrected for extinction, for GRB 990510, at z = 1.619. 
Upper panel: Without subtraction of the host galaxy's 
contribution (the straight line). Lower panel: With the 
host galaxy subtracted, and the CB's AG (the line of 
squares) given by Eqs. (|3l|) and (^). The contribution 
from a 1998bw-like supernova placed at the GRB's red- 
shift, Eq. (0), corrected for extinction, is indicated in both 
panels by a line of crosses. The expected supernova con- 
tribution is too weak to be observed. 



Fig. 14. — Comparisons between our fitted R-band after- 
glow (upper curves) and the observations listed in Table II, 
not corrected for extinction, for GRB 990712, at z = 0.434. 
Upper panel: Without subtraction of the host galaxy's 
contribution (the straight line). Lower panel: With the 
host galaxy subtracted, and the CB's AG (the line of 
squares) given by Eqs. ( plj ) and (p9|). The contribution 
from a 1998bw-like supernova placed at the GRB's red- 
shift, Eq. (0), corrected for extinction, is indicated in both 
panels by a line of crosses. The SN contribution is clearly 
discernible, but a bump at slightly earlier times than that 
of our standard-candle SN1998bw would provide a better 
description. 
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Fig. 15. — Comparisons between our fitted V-band after- 
glow (upper curves) and the observations listed in Table II, 
not corrected for extinction, for GRB 990712, at z = 0.434. 
Upper panel: Without subtraction of the host galaxy's 
contribution (the straight line). Lower panel: With the 
host galaxy subtracted, and the CB's AG (the line of 
squares) given by Eqs. (|3l|) and (^). The contribution 
from a 1998bw-like supernova placed at the GRB's red- 
shift, Eq. (^, corrected for extinction, is indicated in 
both panels by a line of crosses. The SN is clearly dis- 
cernible, but a bump at slightly earlier times than that 
of our standard-candle SN1998bw would provide a better 
description. 



Fig. 16. — Comparisons between our fitted R-band after- 
glow (upper curves) and the observations listed in Table II, 
not corrected for extinction, for GRB 991208, at z = 0.706. 
Upper panel: Without subtraction of the host galaxy's 
contribution (the straight line). Lower panel: With the 
host galaxy subtracted, and the CB's AG (the line of 
squares) given by Eqs. ( pl| ) and (p9|). The contribution 
from a 1998bw-like supernova placed at the GRB's red- 
shift, Eq. (0), corrected for extinction, is indicated in both 
panels by a line of crosses. The SN contribution is clearly 
discernible. 
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Fig. 17. — Comparisons between our fitted R-band after- 
glow (upper curves) and the observations listed in Table II, 
not corrected for extinction, for GRB 991216, at z = 1.020. 
Upper panel: Without subtraction of the host galaxy's 
contribution (the straight line). Lower panel: With the 
host galaxy subtracted, and the CB's AG (the line of 
squares) given by Eqs. (|3l|) and (^). The contribution 
from a 1998bw-like supernova placed at the GRB's red- 
shift, Eq. (0), corrected for extinction, is indicated in both 
panels by a line of crosses. The data show possible evi- 
dence for a SN bump. 



Fig. 18. — Comparisons between our fitted R-band after- 
glow (upper curves) and the observations listed in Table II, 
not corrected for extinction, for GRB 000131, at z = 4.5. 
Upper panel: Without subtraction of the host galaxy's 
contribution (the straight line). Lower panel: With the 
host galaxy subtracted, and the CB's AG (the line of 
squares) given by Eqs. ( plj ) and (p9|). The contribution 
from a 1998bw-like supernova placed at the GRB's red- 
shift, Eq. (0), corrected for extinction, is indicated by a 
line of crosses. The data are too scarce for conclusions to 
be drawn. 
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Fig. 19. — Comparisons between our fitted R-band af- 
terglow (upper curves) and the observations listed in Ta- 
ble II, not corrected for extinction, for GRB 000301c, 
at z = 2.040. Upper panel: Without subtraction of the 
host galaxy's contribution (the straight line) . Lower panel: 
With the host galaxy subtracted, and the CB's AG (the 
line of squares) given by Eqs. ( |3l| ) and (p9|). The contri- 
bution from a 1998bw-like supernova placed at the GRB's 
redshift, Eq. (^, corrected for extinction, is indicated by a 
line of crosses. The expected SN contribution is too weak 
to be observable. 



Fig. 20. — Comparisons between our fitted R-band after- 
glow (upper curves) and the observations listed in Table II, 
not corrected for extinction, for GRB 000418, at z = 1.119. 
Upper panel: Without subtraction of the host galaxy's 
contribution (the straight line). Lower panel: With the 
host galaxy subtracted, and the CB's AG (the line of 
squares) given by Eqs. ( plj ) and (p9|). The contribution 
from a 1998bw-like supernova placed at the GRB's red- 
shift, Eq. (0), corrected for extinction, is indicated in both 
panels by a line of crosses. There is some indication of a 
SN1998bw-like contribution. 
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Fig. 21. — Comparisons between our fitted R-band after- 
glow (upper curves) and the observations listed in Table II, 
not corrected for extinction, for GRB 000926, at z = 2.066. 
Upper panel: Without subtraction of the host galaxy's 
contribution (the straight line). Lower panel: With the 
host galaxy subtracted, and the CB's AG (the line of 
squares) given by Eqs. (|3l|) and (^). The contribution 
from a 1998bw-like supernova placed at the GRB's red- 
shift, Eq. (0), corrected for extinction, is indicated in both 
panels by a line of crosses. A SN1998bw-like contribution 
could not have been seen. 



Fig. 22. — Comparisons between our fitted R-band after- 
glow (upper curves) and the observations listed in Table II, 
not corrected for extinction, for GRB 010222, at z > 1.474. 
Upper panel: Without subtraction of the host galaxy's 
contribution (the straight line). Lower panel: With the 
host galaxy subtracted, and the CB's AG (the line of 
squares) given by Eqs. ( plj ) and (p9|). The contribution 
from a 1998bw-like supernova placed at the GRB's red- 
shift, Eq. (0), corrected for extinction, is indicated in both 
panels by a line of crosses. There are no observations at 
the time a 1998bw-like supernova would have significantly 
contributed. 
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Fig. 23. — Comparisons between our fitted R-band after- 
glow (upper curves) and the observations listed in Table II, 
not corrected for extinction, for GRB 970508, at z = 0.835. 
Upper panel: Without subtraction of the host galaxy's 
contribution (the straight line). Lower panel: With the 
host galaxy subtracted, and the CB's AG (the line of 
squares) given by Eqs. (|3l|) and (^). The contribution 
from a 1998bw-like supernova placed at the GRB's red- 
shift, Eq. (0), corrected for extinction, is indicated in both 
panels by a line of crosses. This is the best fit not including 
the possibility of lensing; and it is terrible. 



Fig. 24. — The same as in Fig. (p3|), but with a den- 
sity profile that jumps from one constant value to an- 
other, at a distance ^ 0.24 kpc from the progenitor. Upper 
panel: Without subtraction of the host galaxy's contribu- 
tion (the straight line). Lower panel: With the host galaxy 
subtracted, and the CB's AG (the line of squares) given 
by Eqs. (|l|) and ^. This time the fit, unhke that of 
Fig. (P3|), is quite satisfactory. A SN1998bw-like contribu- 
tion is necessary. 
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Fig. 26. — The early-time and late-time X-ray AG of 
GRB 970228 in the 2-10 keV range as measured by Costa 
et al. (1997), fitted with Eq. ( [l3| ) for a constant density 
along the CB trajectory. This assumption should be vastly 
inappropriate between ~ 2 x 10"'^ and ~ 0.2 days, during 
which lie and the X-ray AG are expected to diminish by 
two or three orders of magnitude. Coincidentally, there are 
no data there. 
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Fig. 25. — Comparisons between our fitted R-band after- 
glow (upper curves) and the observations, not corrected 
for extinction, for GRB 990123. Upper panel: Without 
subtraction of the host galaxy's contribution (the straight 
line). Lower panel: With the host galaxy subtracted. Un- 
like in the previous ensemble of figures, these AG data 
start from a short time after the GRB. The starting t^^ 
behaviour is that expected if the CBs are moving trough 
a density profile, n c>c r~^, induced by the parent-star's 
pre-SN wind and ejections. 
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Fig. 27. — The early-time and late-time X-ray AG of 
GRB 970828 in the 2-10 keV range as measured by Smith 
et al. (2001) and Yoshida et al. (2001), fitted with Eq. (H) 
for a constant density along the CB trajectory. The same 
comments as in Fig (Eq) apply here. 
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Fig. 28. — The early-time and late-time X-ray AG of 
GRB 971214 in the 2-10 keV range as measured by Dal 
Fiume et al. (2000) fitted with Eq. (||) for a constant 
density along the CB trajectory. The same comments as 
in Fig ( Pq ) apply here. 



Fig. 30. — The early-time and late-time X-ray AG of 
GRB 990123 in the 2-10 keV range (Piro 2000) fitted with 
Eq. ( [43|) for a constant density along the CB trajectory. 
The same comments as in Fig (|2^) apply here. 
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Fig. 29. — The early-time and late-time X-ray AG of 
GRB 980613 in the 2-10 keV range (Pho 2000) fitted with 
Eq. (|^ for a constant density along the GB trajectory. 
The same comments as in Fig (Eq) apply here. 
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Fig. 31. — The early-time and late-time X-ray AG of 
GRB 990510 in the 2-10 keV range as measured by Pian et 
al. (2001) fitted with Eq. (|4^) for a constant density along 
the CB trajectory. The same comments as in Fig (26) 
apply here. 
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Fig. 32. — The early-time and late-time X-ray AG of 
GRB 000926 in the 2-10 keV range as measured by Piro et 
al. (2001), fitted with Eq. (^) for a constant density along 
the CB trajectory. The same comments as in Fig ( p6| ) 
apply here. 
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Fig. 33. — The early-time and late-time X-ray AG of 
GRB 010222 in the 2-10 keV range as measured by In 
't Zand et al. (2001), fitted with Eq. ( |4^ ) for a constant 
density along the CB trajectory. The same comments as 
in Fig (g^) apply here. 
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Fig. 34. — The early-time and late-time X-ray AG of 
GRB 970508 in the 2-10 keV range as measured by Piro 
et al. (1998), fitted with Eq. (|3|) for a - l/r^ plus con- 
stant density along the GB trajectory. The contribution of 
the X-ray line observed before 0.8 d was subtracted from 
the data. This reduces the two observed points between 
0.2 and 1 day by a factor ~ 0.39. The overall result is com- 
patible with an effect that, at late times, is achromatic, 
since the late optical and X-ray AGs are both proportional 
to Uc; see Fig. (gj) for the optical counterpart. 



''■■i \ 

\ : 



10 ■ 10 1 

Time After Burst (days) 



Fig. 35. — The early-time and late-time X-ray AG of 
GRB 991216 in the 2-10 keV range as measured by Corbet 
and Smith (1999) and Takeshima et al. (1999) fitted with 
Eq. ( p3| ) for a density declining like r^^ plus a constant 
density along the CB trajectory. 
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Fig. 36. — Upper panel: A fit to the X-ray afterglow 
of the SN1998bw/GRB 980425 pair. We caU "plateau" 
the slowly-declining late measurements. Lower panel: The 
V-band light curve of the same pair, with the blue "SN" 
curve a fit to the SN by SoUerman et al. (2000), dominated 
after day ^^40 by ^^Co decay. The red "AG" curve is our 
prediction for the CB-induced AG component, as given by 
Eq. (pl|), with the parameters determined from the X-ray 
AG fit in the upper panel. The SN contribution dominates 
up to day ~600. The last point is an HST measurement at 
day 778, that precisely agrees with the (dashed) SN plus 
CB prediction for the total AG. For an earlier version of 
these results, see DD2000a. 
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Fig. 37. — Distribution of a values, from the fits to the 
optical AGs. The prediction is a « 1.1. Binning from a — 
1.06 to 1.26 would have made this distribution look even 
more impressively narrow. The GRBs are labelled by the 
last three digits of their date. 
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Fig. 38. — Distribution of 70 values from the fits to the 
optical AGs. The expectation is 70 ^ 10"^. The GRBs are 
labelled by their last three digits. 
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Fig. 41. — Distribution on E'P^ values, as given by 
Fig. 39.- Distribution of 6 values, in miUiradians, from ^^ Q^ -^ ^^^^ ^^ ^^44 ^^^_ ^he GRBs are labelled by 
the fits to the optical AGs. The GRBs are labelled by their ^^^^^ j^^^ three digits 
last three digits. 
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Fig. 40. — Distribution of Logio[xoo(Mpc)] values, from 
the fits to the optical AGs. The expectation from Eq. ( p8[ ) 
is 0.114. The GRBs are labelled by their last three digits. 



